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As Charles Goodyear discovered in 1839, when he first vulcanised rubber, a macromolecular liquid 
is transformed into a solid when a sufficient density of permanent crosslinks is introduced at random. 
At this continuous equilibrium phase transition, the liquid state, in which all macromolecules are 
delocalised, is transformed into a solid state, in which a nonzero fraction of macromolecules have 
spontaneously become localised. This solid state is a most unusual one: localisation occurs about 
mean positions that are distributed homogeneously and randomly, and to an extent that varies 
randomly from monomer to monomer. Thus, the solid state emerging at the vulcanisation transition 
is an eqmUbrium amorphous solid state: it is properly viewed as a solid state that bears the same 
relationship to the liquid and crystalline states as the spin glass state of certain magnetic systems 
bears to the paramagnetic and ferromagnetic states, in the sense that, like the spin glass state, it is 
diagnosed by a subtle order parameter. 

In this article we give a detailed exposition of a theoretical approach to the physical properties of 
systems of randomly, permanently crosslinked macromolecules. Our primary focus is on the equi- 
librium properties of such systems, especially in the regime of Goodyeax's vulcanisation transition. 
This approach rests firmly on techniques from the statistical mechanics of disordered systems pio- 
neered by Edwards and co-workers in the context of macromolecular systems, and by Edwards and 
Anderson in the context of magnetic systems. We begin with a review of the semi-microscopic formu- 
lation of the statistical mechanics of randomly crosslinked macromolecular systems due to Edwards 
and co-workers, in particular discussing the role of crosslinks as quenched random variables. Then 
we turn to the issue of order parameters, and review a version capable, inter alia, of diagnosing the 
amorphous solid state. To develop some intuition, we examine the order parameter in an idealised 
situation, which subsequently turns out to be surprisingly relevant. Thus, we are motivated to hy- 
pothesise an explicit form for the order parameter in the amorphous solid state that is parametrised 
in terms of two physical quantities: the fraction of localised monomers, and the statistical distribu- 
tion of localisation lengths of localised monomers. Next, we review the symmetry properties of the 
system itself, the liquid state and the amorphous solid state, and discuss connections with scattering 
experiments. Then, we review a representation of the statistical mechanics of randomly crosslinked 
macromolecular systems from which the quenched disorder has been eliminated via an application 
of the replica technique. We transform the statistical mechanics into a field-theoretic representation, 
which exhibits a close connection with the order parameter, and analyse this representation at the 
saddle-point level. This analysis reveals that sufficient crosslinking causes an instability of the liquid 
state, this state giving way to the amorphous solid state. To address the properties of the amor- 
phous solid state itself, we solve the self-consistent equation for the order parameter by adopting 
the hypothesis discussed earlier. Hence, we find that the vulcanisation transition is marked by the 
appearance of a nonzero fraction of localised monomers, which we compute, the dependence of this 
fraction on the crosslink density indicating a connection with random graph theory and percolation. 
We also compute the distribution of localisation lengths that characterises the ordered state, which 
we find to be expressible in terms of a universal scaling function of a single variable, at least in the 
vicinity of the transition. Finally, we analyse the consequences of incorporating a certain specific 
class of correlations associated with the excluded-volume interaction. 

PACS numbers: 61.43, 83. 80. J, 64.70.Pf, 64.60.A , 64.60. C, 36.20, 75.10.N 
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There is probably no other inert substance, 
the properties of which excite in the human 
mind, when first called to examine it, an equal 
amount of curiosity, surprise, and admiration. 
Who can examine, and reflect upon [the properties of] 
gum-elastic, without adoring the wisdom of the Creator! 

Charles Goodyear, Gum-Elastic and its Varieties, 
with a Detailed Account of its Applications and Uses, 
and of the Discovery of Vulcanization (1855) 

I. INTRODUCTION AND OVERVIEW 

In this article we aim to present a theoretical description of the physical properties of systems of macromolecules 
that have been randomly and permanently crosslinked. Our focus will be on the equilibrium properties of such 
systems, especially in the regime of the vulcanisation transition. By the term vulcanisation transition we mean the 
sharp thermodynamic phase transition occurring when the mean density of crosslinks exceeds a certain critical value. 
At this critical crosslink-density, the equilibrium state of the system undergoes a continuous transition: for subcritical 
values the equilibrium state is a liquid state, in which all the macromolecules are delocalised; for supercritical values 
the equilibrium state is an amorphous solid state, in which a nonzero fraction of macromolecules form a macroscopic 
network, and spontaneously become localised, albeit about certain random locations. Thus, our focus will be on the 
spontaneous emergence of the equilibrium amorphous solid state at the vulcanisation transition, and the consequent 
properties of this unusual state. 

It must be emphasised at the outset that the theoretical description of the vulcanisation transition to and properties 
of the amorphous solid state of randomly crosslinked macromolecular networks presented here represents the confluence 
of two pioneering contributions to the theory of condensed matter: the Deam-Edwards theory of a single crosslinked 
macromolecule ^ , and the Edwards- Anderson theory of spin glasses |^ . 

The basic ingredients of the present approach to the physical properties of randomly crosslinked macromolecular 
networks are as follows. We adopt a semi-microscopic description of the macromolecules, in which the detailed 
microscopic chemistry of the macromolecules and solvent (if any there be) feature only to the extent that they 
determine the following effective parameters: the total arclength of each macromolecule, the persistence length (i.e., 
the length of the statistically independent macromolecular segments, which we term monomers), and the excluded- 
volume strength (i.e., the parameter that describes the effective repulsion between monomers). Thus, we regard 
the macromolecules as extended, featureless, flexible linear objects, each capable of exhibiting a large number of 
configurations, and use classical equilibrium statistical mechanics to address the properties of systems composed of 
a thermodynamically large number of such macromolecules. We regard the crosslinks as permanent elements that 
constrain certain randomly chosen monomers to remain adjacent to one another. Thus the crosslinked macromolecular 
system is a system with quenched disorder, in the sense that in addition to the macromolecular freedoms — the so-called 
annealed variables, which undergo equilibrium statistical-mechanical fluctuations — there are additional variables, 
those that specify the detailed realisation of the crosslinking, that do not undergo equilibrium statistical-mechanical 
fluctuations. Instead, these variables — the so-called quenched random variables — vary only between realisations of the 
physical system. We treat these quenched random variables statistically, too, but account for their quenched nature 
by invoking the replica technique. 

We now give an overview of our basic strategy for determining the physical properties of randomly crosslinked 
macromolecular networks. 

We characterise the plausible equilibrium states of the system — liquid, globule, crystalline solid, amorphous solid — 
in terms of an order parameter, designed to discriminate between these states. This order parameter is a vastly 
more intricate object than the order parameters that arise, say, in the study of ferromagnetism or even spin glasses, 
and one must explore enormously larger spaces to find its equilibrium value. By investigating a simple caricature 
of the amorphous solid state we are, however, able to identify a physically well-motivated scheme for parametrising 
the amorphous solid state order parameter at a manageable level: via a single number — the fraction of spatially 
localised monomers — and a normalised probability distribution — the statistical distribution of localisation lengths of 
the localised monomers. 

We focus on the free energy and the order parameter for the system of interacting macromolecular freedoms subject 
to the crosslinking constraints. Application of the replica technique to these quantities allows for the elimination of 
the quenched random variables; the price for this elimination is the introduction of an unusual and rather compli- 
cated effective coupling amongst the replicated macromolecular freedoms. Our scheme for parametrising the order 
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parameter leaves intact the permutation symmetry amongst the repHcas. Next, we introduce a certain stochastic 
field, the argument of which is replicated (real or wave-vector) space, via which we transform the replicated macro- 
molecular description into a field-theoretic one. In this representation, the individual macromolecules are coupled to 
one another only indirectly, via their coupling to the fluctuations of the stochastic field, although the replicas of any 
given macromolecule remain directly coupled to each other. Via the introduction of a suitable external potential we 
demonstrate the connection between the order parameter and the stochastic field. 

In order to elucidate the physical properties of the system we adopt a mean-field approach, which amounts to our 
approximating, via the saddle-point method, the averages over the stochastic field in the field-theoretic expressions for 
the free energy and order parameter. The state of the physical system then follows from the form of the appropriate 
stationary value of the stochastic field or, equivalently, from the form of the self-consistent value of the order parameter. 
For mean crosslink densities smaller than a certain critical value, of order one crosslink per macromolecule, there is 
only one stationary value, which is elementary, and the corresponding state is the liquid state. For supercritical 
crosslink densities the appropriate stationary value is most definitely not elementary, and the corresponding state is 
the amorphous solid state. At the critical crosslink density the system undergoes a continuous vulcanisation transition 
from the liquid state to the amorphous solid state. 

To determine the properties of the amorphous solid state itself, we hypothesise that the self-consistent value of the 
order parameter lies within the family of order-parameter values reachable via our (severely restrictive but nevertheless 
physically plausible) parametrisation. Quite remarkably, this is indeed the case: our parametrisation does not merely 
yield a variational approximation to the amorphous solid state. Instead, although one has no a priori reason to 
suppose that it should, it permits an exact mean-field description of randomly crosslinked macromolecular networks to 
be constructed. What emerges is an amorphous solid state characterised by a nonzero fraction of localised monomers. 
The precise value of this fraction depends on the crosslink density, and vanishes continuously at the transition and in 
the liquid state. This fraction depends on the crosslink density in a manner identical to that found in random graph 
theory and percolation. The state is further characterised by a crosslink-density-dependent distribution of localisation 
lengths, which quantifies the manner in which the localised monomers have become localised around their random mean 
positions. The typical localisation length diverges continuously at the transition and in the liquid state. In the vicinity 
of the transition, the distribution of localisation lengths has a scaling form governed by a universal function, which we 
compute. To date, we have been unable to obtain conclusive results for the distribution of localisation lengths far from 
the amorphous solidification transition. The reason for this is purely technical: at a certain stage in the development 
we employ a perturbative calculation, in which the small parameter is the characteristic inverse localisation length, 
measured in units of the radius of gyration of a single, noninteracting macromolecule, this parameter being zero in 
the liquid state, and small in the amorphous solid state only in the vicinity of the transition. 

In addition to developing a purely mean-field description of the transition to and properties of the amorphous solid 
state, we incorporate a class of correlations associated with the excluded-volume interactions. We accomplish this by 
treating a certain sector of the stochastic field at the gaussian level, rather than the saddle-point level. What emerges 
from this treatment is merely a particularly simple finite renormalisation of a specific parameter of the mean-field 
theory, at least in three spatial dimensions. 

Approaches based on percolation theory |^ are sometimes pursued in the context of the modelling of randomly 
crosslinked macromolecular networks; see Ref. However such approaches entail only a single statistical ensemble, 
and hence cannot treat the equilibration of thermal (i.e., annealed) freedoms in the presence of quenched freedoms 
(crosslinks, in the present case). If the percolation ensemble is interpreted as describing the quenched freedoms 
then it can account for the appearance of an infinite network at a critical crosslink density, but it cannot account 
for the thermal fluctuations, which determine the physical properties of the liquid and amorphous solid states. In 
particular the degree of localisation in the amorphous solid state and the entropic elasticity of the network cannot be 
calculated from first principles in a percolation-based model, which does not distinguish between thermal and quenched 
freedoms. If, on the other hand, the percolation ensemble is interpreted as describing the thermal fluctuations of the 
macromolecular freedoms then it may serve as a model for chemical gelation, in which crosslinks continuously form 
and break up, but it cannot account for permanent, quenched crosslinks. 

An approach to the properties of randomly permanently crosslinked macromolecular systems that accounts for 
both thermal fluctuations and quenched disorder has been introduced by Edwards and co-workers [Qj^J^ . In addition 
to formulating the statistical mechanics of a single randomly crosslinked macromolecule using the replica technique, 
Deam and Edwards investigated the elastic properties of the system by invoking a variational method Ball and 
Edwards investigated the incorporation of correlations, and Ball also studied certain aspects of many- macromolecule 
systems 0. 

The point of view that the solid state of randomly crosslinked macromolecular networks represent an unusual, 
equilibrium amorphous solid state of matter was introduced and developed in detail in Refs. In addition, 

the order parameter appropriate for detecting this amorphous solid state was introduced and explained in these 
papers, and the instability of the liquid state with respect to the formation of an amorphous solid, induced by 
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sufficient crosslinking, was also identified there. Extending these ideas, the transition to and certain properties of the 
amorphous sohd state, including its elastic properties, were obtained within the context of a variational mean-field 
approach in Ref. This variational approach was improved through an exact mean- field approach, the central 

elements of which were reported in Ref. [[l3[ . The present article is intended to provide a fairly complete exposition 
of this approach to the statistical mechanics of randomly crosslinked macromolecular networks, including a detailed 
account of the work reported in Refs. ||l^ , [l3|] . 

By adopting the Deam-Edwards approach ^ and the amorphous solid order parameter picture of Refs. [^,^, 
Panyukov has made a number of contributions to the theory of well-crosslinked macromolecular networks To 
accomplish these, he considers a single linear macromolecule, randomly crosslinked, which is intended to represent a 
physical network of many randomly crosslinked macromolecules. These contributions, which are reviewed in Ref. [p^ , 
are based on the introduction of a well-known free-field representation of a random walk p6| , the resulting field 
being treated at the mean-field level, along with an additional replica limit. Focusing, as it does, on the statistical 
mechanics of a single, well-crosslinked macromolecule, Panyukov's approach is unsuitable for developing a theory of 
the transition between the liquid and amorphous solid states, which occurs at crosslink densities of order one per 
macromolecule . 

The present approach to the physics of randomly crosslinked macromolecular systems possesses the following dis- 
tinctive virtues. First, both thermal freedoms (i.e., the macromolecular positions) and quenched freedoms (i.e., the 
crosslink locations) are incorporated, and handled appropriately, in contrast with percolative pictures. The replica 
technique provides the tool for accomplishing this. In fact, the percolative picture emerges from the present approach 
in the form of statistical information concerning the presence of localised macromolecules. However, the present 
approach is considerably richer, additionally yielding statistical-mechanical information about the (thermally fluc- 
tuating) macromolecular system. In particular, it allows a unified treatment of liquid and amorphous solid states. 
Second, the physical many-macromolecule character of the system is maintained, in contrast with approaches that 
consider instead the properties of a single macromolecule. Especially in the vicinity of the solidification transition, 
where the number of physical crosslinks is of order one per macromolecule, this is particularly significant, and al- 
lows us to develop a theory of the solidification transition. Third, the present approach leads directly to an order 
parameter for the amorphous solid state, which is related to that of spin glass physics. The order parameter has a 
natural, physical interpretation, which facilitates the hypothesising of an appropriate form for it. Fourth, the physical 
freedoms, viz., the macromolecular configurations, appear directly throughout the development, not being exchanged 
for some problem-specific formal representation. We can, therefore, be confident that the macromolecular character 
of the system is retained, especially when approximations are made. Indeed, the entire approach is very robust, so 
that, in addition to being of interest in the context of vulcanised macromolecular systems, it can readily be extended 
to address a wide range of other physical systems, such as crosslinked manifolds cndlinked systems of flexible, 
semi- flexible and rigid macromolecules , continuous random network models of structural glasses [p^ , and proteins. 
In addition, it should prove possible to extend the present approach to address issues of dynamics. Moreover, looking 
beyond the algebraic details, one sees a theoretical superstructure that is rather natural, direct and perhaps even 
conventional, at least from the point of view of statistical fleld theory. Fifth, the use of the present approach has 
primarily been restricted to the mean-field level of approximation. However, it is not a variational mean-field theory: 
instead the relevant saddle-point is determined exactly. This makes the approach a particularly promising starting 
point for future developments, such as the elucidation of the elastic properties and of the role of thermal fiuctuations. 
It should be mentioned that the full stability of the saddle-point that we have determined has not yet been established. 
One current shortcoming of the present approach arises from the technical difficulty of computing properties in the 
regime of high crosslinking (e.g., deep in the amorphous solid state). This does not limit the scope of our primary 
aim, viz., to explore the vicinity of the transition to the amorphous solid state. A second shortcoming is the inability 
of the approach, at least in its present formulation, to respect the interlocking of closed loops of macromolecules 
(discussed in Sec. HE) that crosslinking can induce. 

The present article is organised as follows. In the present, introductory, section we provide an overview of the article. 
In Scc.g we discuss the basic elements of the model of macromolecular systems on which the present approach is based, 
including the level of description of macromolecular configurations, the Edwards measure for their statistical weights, 
and the notion of crosslinks as quenched random variables. We also discuss the partition function, free energy and 
issues of indistinguishability, along with the statistical characterisation of the crosslinks, and the notion of disorder 
averages of certain physically relevant quantities. In Sec. Ill we develop the general subject of order parameters 
appropriate for the amorphous solid and other states, discussing the properties that such order parameters should 
possess. We explore a simple scenario for the amorphous solid state, which provides physical motivation for a certain 
specific hypothesis that we make for the form taken by the amorphous solid order parameter in the amorphous solid 
state. At this stage we introduce the concept of gel and sol fractions and the statistical distribution of localisation 
lengths associated with localised monomers. We also analyse the symmetry properties of the ordered state, and 
discuss the connection between the order parameter and the elastic scattering of neutrons. In Sec. IV we address 
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the statistical mechanics of randomly crosslinked macromolecular networks, invoking the replica technique in order 
to eliminate the (quenched random) crosslink variables. We also introduce a suitable replica-Helmholtz free energy, 
which is dependent on a convenient generalised external potential. In Sec. ^ we reformulate the statistical mechanics 
of randomly crosslinked macromolecular networks in field-theoretic terms by introducing a certain stochastic field, 



which is closely related to the amorphous solid state order parameter. In Sec. VI we explore the properties of the 
resulting field theory within the context of a natural mean-field approximation. We exhibit the instability of the 
liquid state, and we compute the free energy and self-consistent order parameter in the vicinity of the transition. 



We also describe the characteristics of the amorphous solid state that emerge from this approach. In Sec. VII 



incorporate certain correlations associated with density-sector fluctuations, and demonstrate that the resul ts of the 
previous sections of this article are robust with respect to the incorporation of these fluctuations. In Sec. VIII we make 
some concluding remarks. We have organised this article so that the main text is, to a large degree, unencumbered 
by lengthy mathematical details. Wherever possible such details have been relegated to one of eleven appendices, in 
which we demonstrate in full detail how the results of the main text are established. 

Perhaps the most significant property to emerge at the vulcanisation transition is rigidity with respect to shear 
deformations, i.e., elasticity. In a forthcoming article ||20|| we shall report in detail on the extension of the present work 
to the investigation of the elastic properties of the amorphous solid state of randomly crosslinked macromolecular 
networks. 



II. MODEL OF THE MACROMOLECULAR SYSTEM 

A. Macromolecular system prior to crosslinking 

We consider a system consisting of a large number N of long, flexible macromolecules, initially identical, and 
contained in a large d-dimensional hypercube of volume V. The macromolecules are characterised by their com- 
mon arclength L and (weakly temperature-dependent) persistence length ^(<C L), so that the number of effectively 
statistically independent segments comprising each macromolecule is of order L/£ 3> 1. Semi- microscopic spatial con- 
figurations of the system are characterised by the collection of spatial configurations of the macromolecules {'Ri{a)}fLi, 
in which Ri(a') is the d-dimensional position vector of the monomer an arclength distance a from a specific end of 
macromolecule i, the (discrete) macromolecule index i ranging from 1 to iV and the (continuous) arclength variable 
a ranging from to L. 

It is convenient to exchange the dimcnsionful position vector R and arclength a for dimensionless versions c and s 
via the transformation 



R,(ct) = ^/£L/d c,;(s), (2.1a) 
a = Ls. (2.1b) 

1 /2 

Thus, we shall be measuring spatial distances in units of (^£L/d) {i.e., the root mean squared end-to-end distance 

of a free macromolecule divided by Vd), and arclength distances in units of the total arclength L. We shall measure 
energies in units such that k^T is unity. 

At the level of the present semi-microscopic description, and prior to the incorporation of the effects of either 
monomer-monomer interactions or crosslinks, we account for the connectivity of the constituent macromolecules by 
adopting the Wiener measure pl|-p3[ , in terms of which the statistical weight of the configuration of the system 
{ci{s)}^i is proportional to exp ( — H^), where 



i=i 



W.I ^ 
ds 



(2.2) 



The subscript 1 on anticipates the introduction of replicas of the system, which we shall need to make below. 
We shall often need to consider normalised expectation values taken with respect to the Wiener measure, which we 
shall denote by the angle-bracket pair (• • defined by 



J Vc exp (-Fi^) 
^ f Vc exp i-H^ 



(2.3) 
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where the dots represent an arbitrary function of the configuration of the system, and the measure 



N 



I?c = J| Jl dc,{s) 



(2.4) 



1=1 0<s<l 



indicates functional integration over aU spatial configurations of the system, i.e., over all configurations of the N 
macromolecules. The subscript 1 on (• • •) also indicates that the average is taken only over the configurations of a 
single copy of the system, also anticipating the introduction of replicas. 

We account for monomer- monomer interactions in a phenomenological manner, by augmenting the Wiener measure 
with an additional factor that has the effect of suppressing the statistical weight of configurations in whi ch p airs 
of monomers occupy common regions of space [pi[-23|. To this end, we replace the Wiener measure, Eq. (2.2), by 
the Edwards measure, in terms of which the statistical weight of the configuration {ci{s)}^i is proportional to 
exp ( — Hf^ , where 



1 ^ 



ds 



d 2 



.2 N 



ds 



ds'(5('^)(c,(s)-c,,(s'))- 



(2.5) 



Here, (5^'^^(c) is the d-dimensional Dirac (5-function, and the dimensionless (real, positive) parameter charac- 
terises the strength of the suppression of statistical weight due to the (repulsive) excluded- volume interaction between 
monomers pl|-p3|. The excluded- volume interaction is suitably modified so as to exclude interactions between ad- 
jacent monomers on a common macromolecule (i.e., monomers for which |s — s'| < i/L). The expression for the 
Edwards measure in terms of dimensionful variables is given in footnote p5|] . The system can be regarded as a melt 
of macromolecules, in which case the interaction parameter is intended to account for the monomer- monomer 
interaction. Alternatively, it can be regarded as a solution of macromolecules dissolved in a good solvent, in which 
case A^ is intended to represent the effective monomer-monomer interaction (i.e., the bare interaction renormalised 
by the monomer-solvent and solvent-solvent interactions, the solvent degrees of freedom having been integrated out). 
In both cases, A^ is weakly temperature-dependent. Even at the level of mean-field theory, the excluded- volume inter- 
action plays a crucial role: it partially compensates the effective monomer- monomer attraction due to the crosslinks 
in just such a fashion as to maintain the macroscopic homogeneity of the system whilst allowing for the possibility of 
transition from the liquid to the amorphous solid state. 

We shall need to consider normalised expectation values taken with respect to the Edwards measure, which we shall 
denote by the angle-bracket pair {• ■ ■)f, defined by 



J Vc exp (-iff) • ■ 
J Vc exp (-iff) 



(2.6) 



where the dots represent an arbitrary function of the configuration of the system, and once again Vc indicates 
functional integration over all configurations of the system. Again, the subscripts 1 on Hf and (• • anticipate the 
introduction of replicas. 

It should be noted that neither the Wiener measure nor the Edwards measure explicitly breaks translational or ro- 
tational symmetry: the statistical weight of a configuration remains unchanged if all the monomers are simultaneously 
translated through a common amount or rotated through a common angle about a common axis. 



B. Crosslinks as quenched random variables 



Our aim is to address the statistical mechanics of thermodynamically large systems of macromolecules into which a 
large number of crosslinks have been permanently introduced at random. Each crosslink has the effect of constraining 
two randomly selected monomers, the locations of which were kinematically independent prior to the introduction of 
the crosslink, to occupy a common spatial location. Thus the effect of the crosslinks is to eliminate from the ensemble 
of configurations of the system all configurations that do not obey the entire set of random constraints enforced by 
the crosslinks. Our task is therefore to address the statistical mechanics of macromolecular systems in the presence 
of a large number of random constraints. 

A specific realisation of the crosslinking is fully described by specifying which randomly selected pairs of monomers 
are connected by each crosslink, i.e., that the crosslink labelled by the index e serves to connect the monomer at 
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arclength Se on macromolecule ie to the monomer at arclength s'^ on macromolecule i'^, for e = 1, . . . , M, with M 
being the total number of crosshnks. Thus, only those configurations that satisfy the constraints 

(se) = (4), (with e = 1, . . . , M) (2.7) 

are retained in the ensemble. It should be noted that these constraints do not explicitly break translational symmetry. 

In principle, of course, neither the crosslinks nor the integrity of the macromolecules are truly permanent. However, 
in many physical realisations of crosslinked macromolecular systems there is a very wide separation between the time- 
scale required for the crosslink-constrained macromolecular system to relax to a state of thermodynamic equilibrium 
and the much longer time-scale required for either the crosslinks or the macromolecules to break. For such systems, 
and it is such systems that we have in mind, the crosslinks and the macromolecules should be regarded as permanent, 
so that the number and identity of the monomers participating in crosslinks, {ie, Se;i'^, s'^J^L^, should be treated 
as nonequilibrating {i.e., quenched) random variables. The unconstrained macromolecular freedoms are regarded as 
reaching equilibrium in the presence of fixed values of the quenched variables. Thus, it is a meaningful task to address 
the equilibrium statistical mechanics of permanently crosslinked macromolecular systems. 

It should be remarked that the relative statistical weights of the configurations that do satisfy the crosslinking 
constraints are hypothesised, at least a priori, to be unaffected by the introduction of crosslinks. That is, the statistical 
weights are proportional to exp (— _ff;p) for configurations satisfying the crosslinks and zero otherwise. However, as 
we shall see in detail below, for a sufficiently large density of crosslinks the translational and rotational symmetry of 
the equilibrium state of the system is spontaneously broken. That is, in a given (pure) state the statistical weights 
of configurations that are translations and rotations of one another are no longer identical, and thus localisation 
can arise. Indeed, only one member of a family of translated and rotated configurations has a nonzero weight in a 
given (pure) state. The associated transition to an amorphous solid state is precisely the transition on which we are 
focusing. We remark that in the present context of amorphous solidification, translational and rotational symmetry 
are spontaneously broken in an unusual sense, in that they remain fully intact at the macroscopic level. 

A second mechanism that leads to the violation of the hypothesis mentioned in the previous paragraph arises 
because sufficient crosslinking is liable to give a topological character to the system of macromolecules, at least in 
three spatial dimensions, in the sense that for a given set of crosslinks there will be families of configurations allowed 
by the crosslinks that are mutually inaccessible. We mean by this that, because of the possibility of interlocking closed 
loops formed by macromolecules, there will be families of configurations between which the system cannot continuously 
evolve without the necessity either of the breaking of at least one crosslink or the passage of one monomer through 
another. We distinguish between constraints arising indirectly from crosslinking via the interlocking of closed loops 
and constraints arising directly from the crosslinks themselves by referring to the former as anholonomic constraints 
and the latter as holonomic constraints. In principle, a statistical-mechanical approach should incorporate, at most, 
those configurations that are mutually accessible, i.e., should respect both holonomic and anholonomic constraints. 
The theory presented here treats the holonomic constraints as quenched but the anholonomic constraints as annealed, 
therefore not incorporating the latter. We know of no explicit semi-microscopic strategy that is capable of handling 
the anholonomic constraints. 



C. Partition function 



We define the statistical-mechanical partition function Z{{ie, Se', i'e, ^eleii) that characterises the crosslinked system 



M 



M g 



(2.8) 



The pro duct of Dirac (5- functions serves to remove from the s um over configurations implicit in the angle brack- 
ets (2^) any configuration that fails to satisfy the constraints (2/7) |2^, the remaining configurations contributing 
with weights given by the Edwards measure (2.5). With this definition, Z{{ie,Se',i'e,s'f,}^i) is normalised relative 
to the uncrosslinked system. Consequently, the free energy derived from this partition function will in fact be the 
increase in free energy that arises upon crosslinking. This allows us to focus on the implications of crosslinking rather 
than the properties of the uncrosslinked system. 

At first sight, the quantity Z({ie, Se; ig, Sg}^]^) in Eq. (2.S), which we are calling the partition function, appears 



to be precisely the physical partition function of the crosslinked system, at least relative to that of the uncrosslinked 
system. Howev er, fo r a straightforward reason associated with the notion of indistinguishability, a reason that we 
discuss in Sees. [ID and HF, Z({ie, Sg; ig, Sg}^i) as defined in Eq. (^) is not quite the correct definition of the 
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physical partition function. However, as we shall see, the partition function Z({ze, Se! ig, Seleii) ^i^^ turn out to be 
adequate for our purposes. 

D. Indistinguishability 

As first pointed out by Gibbs the (configurational aspect of the) physical partition function for systems 

involving one or more species of identical constituents is to be found by summing over all configurations of the system 
whilst ignoring the issue of the distinguishability of the constituents, and subsequently dividing by an appropriate 
combinatorial factor to account for the indistinguishability of the constituents. This strategy compensates for the 
over-counting of configurations that has arisen from the neglect of indistinguishability. 

What are the implications of indistinguishability in the present context? For the case of the system of A'' identical 
uncrosslinked macromolecules, the appropriate factor is Nl, and thus the physical partition function is given by 

^ exp {-Hf) . (2.9) 

If, for the case of the crosslinked system, the appropriate factor were also Nl (which it is not) then its physical 
partition function would be given by 

M 

— J Vce^p{~Hf)l[6^'\cM-c,yj), (2.10) 



and thus the relative physical partition function would indeed be given by Eq. (2.8), the factors of 1/A^! in the 
numerator and denominator cancelling. 

However, the process of crosslinking alters the system from one that comprises N copies of a single species of 
identical elements. Instead, the crosslinked system will contain a variety of species, such as macromolecules that do 
not participate in any crosslinks, as well as clusters of macromolecules of many types. By clusters we mean assemblages 
of macromolecules that are (directly or indirectly) connected by crosslinks or interlockings, and therefore cannot be 
separated by arbitrary distances. Examples of clusters include pairs of macromolecules that participate in a single 
crosslink, that crosslink being located between some specific pair of arclength locations (say (s,s') = (0.12,0.57)), 
single macromolecules crosslinked to themselves at some specific pair of arclength locations, triplets of macromolecules 
connected by two specifically located crosslinks, as well as more complicated species such as pairs of self-crosslinked 
macromolecules interlocking one another. 

Let us label the various possible cluster species by the index a — 1,2,3,..., and let a = label the uncrosslinked 
macromolecule species. Then, for a specific realisation of the disorder {i.e., the crosslinks and the interlockings) let 
the number of uncrosslinked macromolecules be vq, and the number of clusters of species a be I'a [^8| . Then the 
incorrect combinatorial factor of Nl should be replaced by the correct factor of Oa^a'' ^^^^ factor varying across 
disorder realisations. The physical partition function for a given realisation of the system is then given by 

/Pcexp(-i7f)n'^^'Hc«e(se)-c,,«)). (2.11) 

This correction of the combinatorial factor is mirrored by the absence, due to the constraints, in the summation over 
system configurations of those configuration in which macromolecules participating in a cluster are widely separated, 
which results in the loss of volume factors. Together, the corrected combinatorial factor and the loss of volume factors 
conspire to yield a thermodynamic free energy that is properly extensive. 

In common with much work on the physics of disordered systems, we shall not focus on the statistical mechanics 
of a system having a particular realisation of the disorder. Instead we shall take a probabilistic approach, focusing 
on the typical properties of randomly crosslinked macromolecular systems. To do this, we shall need to consider the 
statistical distribution of crosslink locations. In fact, we shall also allow the number of crosslinks to fluctuate across 
realisations. 

E. Deam-Edwards crosslink distribution 

To compute physical quantities characterising the system of randomly crosslinked macromolecules for a specific 
realisation of the large set of quenched random variables {ie, Se', i'f,, Sg}*!]^ is, of course, neither possible nor particularly 
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useful. Instead we shall focus on typical values of physical quantities, constructed by suitably averaging them over 
the quenched random variables. To perform this averaging we shall need to choose a probability distribution that 
assigns a sensible statistical weight PM{{ie, Se;i'e, s'eleii) to each possible realisation of the number M and location 
{?e, Se; Seleli ^hc crossliuks. FoUowing an elegant strategy due to Deam and Edwards we assume that the 
normalised crosslink distribution is given by 

M!(exp(^ E / dsj^ d.'5('^)(c,(s)-c,(s'))))^ 



where Z(^{ie,Se',i'eTs'^}eLi) is given by Eq. (2.8), and can be regarded as probing the equilibrium correlations of the 
underlyinguncrosslinked liquid |^9|] . Such correlations were omitted from the crosslink distribution in certain previous 
works |&-pp6|, which led to difficulties in obtaining a quantitative description of the amorphous solid state. It is 
not, at present, clear whether this omission is significant for the liquid state. 

The Deam-Edwards distribution can be envisaged as arising from a realistic vulcanisation process, in which crosslinks 



are introduced simultaneously and instantaneously into the liquid state in equilibrium 1 30 1 . Specifically, it incorporates 
the notion that all pairs of monomers that happen (at some particular instant) to be nearby are, with a certain 
probability controlled by the crosslink density parameter js^, crosslinked. Thus, the correlations of the crosslink 
distribution reflect the correlations of the uncrosslinked liquid, and it follows that realisations of crosslinks only 
acquire an appreciable statistical weight if they are compatible with some reasonably probable configuration of the 
uncrosslinked liquid. This good feature of the Deam-Edwards distribution is compatible with the random, space-filling, 
frozen liquid, nature of the equilibrium amorphous state that is achieved upon sufficient crosslinking. 

We allow the number of crosslinks to fluctuate in a quasi-Poisson manner, controlled by the parameter fjl^ . All that 
we shall need to know about is that the mean number of crosslinks per macromolecule, which we denote by [M]/N, 
is a smooth, monotonically-increasing function of that can, in principle, be determined using the distribution V 
p^. We remark that the control parameter /i^ appears in Eq. (2.12) divided by N/V. This factor is simply the 
(dimensionless) density of macromolecules, which is an intensive quantity. As we shall see, this choice leads to an 
equation of state that does not depend on the density of macromolecules, at least at the level of mean-field theory. 
We also remark that no delicate scaling of the control parameter is needed to achieve a good thermodynamic limit, 
in contrast with the case of the Sherrington-Kirkpatrick spin-glass model [p2| . 



As discussed in Sec. [IB, at least in three dimensions crosslinking confers anholonomic topological constraints on the 
network as well as holonomic ones. Thus, the statistical-mechanical tool for constructing the crosslink distribution 
is not entirely correct. In principle, crosslink-realisations should be labelled not only by {ie,Se',i'e,s'^}eLi, j-e., by 
the number and arclength locations of the crosslinks, but also by the precise topology of the realisation, i.e., by 
the manner in which the macromolecules thread through the closed loops made by one another. Then the statistical 
weight attributed to a crosslink-and-topology realisation would be better modelled as arising from those configurations 
of the underlying eq uilibr ium liquid that not only satisfy the holonomic constraints but also the anholonomic ones. 



As remarked in Sec. II B, no mathematical tool yet exists for accomplishing this refinement analytically. In other 
words, we are treating the random topology of the system as annealed rather than quenched. 

One should pause to notice the striking feature that at the heart of the Deam-Edwards crosslink distribution is the 
partition function Z(^{ie, Se', i'e, Seleli) of the crosslinked system, i.e., the crosslink distribution is itself proportional to 
the partition function, the logarithm of which it is to be used to average. This fact gives the development a structure 



that is rather appealing, at least from the point of view of form. This will become especially apparent in Sec. IV in 
the context of the replica technique, in which this distribution is generated via an additional {i.e., zeroth) replica, 
the permutation aspect of the symmetry of the theory thereby being enlarged from the permutation group 5„ to 
Sn+i, where n is the number of replicas p3| , ^ . There is, however, no physical basis for restricting attention solely to 
crosslink distributions generated by the partition function identical to that of the crosslinked system. For example, one 
might imagine crosslinking at a different temperature or solvent quality, which would break the symmetry between 
the crosslink distribution and the partition function of the crosslinked system; then, in the context of the replica 
technique, the permutation aspect of the symmetry of the theory would remain iS„. 



F. Disorder averages and symmetry factors 

How are we to use the Deam-Edwards crosslink distribution? As is well known, it is generally inappropriate in 
disordered systems to average the partition function itself over the quenched random variables, as this would amount 
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to treating the quenched random variables as annealed variables {i.e., equilibrated variables having the same status 
as the variables describing the configurations of the system that can be accessed during equilibrium fluctuations). 
Rather, it is thermodynamically extensive or intensive quantities, such as the free energy or the order parameter, that 
should be averaged over the quenched random variables [^5| . To illustrate this point, consider the free energy relative 
to that of the uncrosslinked system, — In Z. (Recall that we are measuring energies in units such that fceT = 1.) Then 
the disordered-average of the free energy per macromolecule (relative to the free energy of the uncrosslinked system) 
per space dimension, which we denote /, is given by 



dNf 



\nZ{{i 



(2.13) 



where the square brackets indicate a disorder-average, viz., 



oo „i „1 
[OM{{le,Se;i'„s'jtil)]=Wo+Yl / dsi---dsM ds[ ■ ■ ■ ds'j^ 

N N N N 

- ^ ^ VM{{ie,Se;ie:S'J^Ll)OM{{ie,Se;i'e:S'jtil), 



(2.14) 



il=l iA/ = li', =1 



where OM{{ie, Se;i'e, s'el^i) is an arbitrary function of the reahsation of crosshnks. The average over the locations 
of the crosslinks excludes realisations of the disorder in which two positions on the same macromolecule located closer 
than a persistence length participate in crosslinks. This can be accomplished by suitably cutting off the arclength 
integrations. Inasmuch as Z is not stri ctly speaking the physical partition function, / is not strictly speaking the 
physical free energy. In fact, from Eq. (2.11) we know that the physical partition function (normalised with respect 
to that of the uncrosslinked system) is given by 



Z{{te,Se;i',,s'jtU 



Nl 



Z{{ie, ' 



Thus, for the disorder-averaged physical free energy /, which is given by 

-dNf^[lnZ{{t,,Se;i:,s'j',U)] , 

we obtain 

- dNf = -dNf + In N\ - \^\nY[val 

a 

Now, / is an intensive quantity. However, A/, defined by 



J J .1 



In 



Nl 



(2.15) 

(2.16) 
(2.17) 

(2.18) 



is in general proportional to In (for large N). The constant of proportionality is, in general, difficult to compute: it 
will, however, be a small number for the case of lightly crosslinked systems, increasing to d~^ in the high-crosslinking 
limit, for which all macromolecules are connected to a single cluster. (For the uncrosslinked system A/ = 0.) Thus, / 
(= / — A/) contains a term proportional to the logarithm of the size of the system, i.e., is not intensive. Despite this 
unphysical feature of /, it is / that we shall be computing, rather than /, because our inability to compute [ In Yla ^a'] 
precludes us from computing /. However, the physical properties of the system, such as the order parameter, e.g., 
are determined by certain disorder-averaged quantities that we shall show to be insensitive to the indistinguishability 
factor Y[a ^a-, and which can thus be computed in the present approach. 



III. ORDER PARAMETER FOR THE AMORPHOUS SOLID STATE 



A. General properties of the order parameter 

We now discuss a certain order parameter constructed with the intention of distinguishing between equilibrium 
states that are liquid (in which the monomers are all delocalised), crystalline solid (in which a nonzero fraction are 
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localised in a spatially periodic fashion) , globular (in which the monomers have condensed within a spatial subvolume 
of the system), and amorphous solid (in which a nonzero fraction are localised in a spatially random fashion) 
This order parameter is an extension, mutatis mutandis, of the order parameter introduced by Edwards and Anderson 
in the context of a class of amorphous magnetic systems known as spin glasses . 

For a specific realisation of the crosslinks {i.e., prior to disorder averaging), the appropriate order parameter is 
given by 



Try] / ds(exp(iki •Ci(s)));^(exp(ik2 
^=l -^0 



(exp [i 



5))>, 



(3.1) 



for (7 = 1,2,3,..., none of the d-dimensional wave vectors {k^, . . . , k^ } being zero. The angle brackets (• • indicate 
an average over the equilibrium state in question for a particular realisation of the disorder, indicated by the subscript 
X. Such equilibrium states may correspond to situations in which the translational symmetry of the system is 
spontaneously broken, in which case they are not ergodic. However, we shall not dwell here on the possibility of 
further ergodicity-breaking (e.g., of the type commonly associated with the concept of replica-symmetry breaking; 
see Ref. [Q). This restriction is consistent with the results presented below. For a discussion of ergodicity-breaking 
in systems of crosslinked macromolecular networks, see Refs. |Hll|,||. 
The disorder-averaged order parameter is denoted by 



N ^1 



Tr'^ ds (exp (ik^ ■ c,(s)))^(exp (ik^ ■ c,{s)))^ ■ ■ ■ (exp {{k^ ■ Ci(s)))^ 



(3.2) 



For any particular positive integer g, this order parameter may be regarded as the g'^ moment of the distribution of 
random static density fluctuations A/'({pk}) (see Ref. Q), which is defined by 



A/'({Pk}) ^ 



'5c(pk - (exp (ik- Ci(s)))x) 



(3.3) 



where Y[k denotes the product over all d-vectors k in the half-space given by the condition k • n > for a suitable 
unit d- vector n, and the Dirac (5-function of complex argument Sc{z) is defined by Sc{z) = S{Kez) SiJraz), where Rez 
and Im2; respectively denote the real and imaginary parts of the complex number z. Thus, 

/ JJ^dRepkrflmpk J^{{pk}) Pk^Pk-^ ■ ■ ■ Pks 

T7 / (exp (ik^ ■ Ci(s)))x(exp (ik^ ■ c,(s)))x • • • (exp {ik^ ■ c,(s))), 
^=l -^0 



(3.4) 



To see why formula ( |3.lD is indeed an order parameter appropriate for distinguishing between liquid, crystalline, 
globular and amorphous solid states, let us examine its qualitative properties. First, suppose that the state is liquid. 
Then each monomer (z, s) is to be found, with equal probability, in the vicinity of any location in the container. 
Consequently, the equilibrium expectation value of its density ((^('^^(r — Ci{s)))^ is the constant V^^, and the Fourier 
transform ( exp (ik ■ Ci{s)))y^ vanishes (except for the trivial case of k = 0). Thus, for a liquid state the order 
parameter ( |3.l| ) vanishes, all terms in the summation over monomers vanishing identically. This corresponds to a 
state having full translational and rotational symmetry. 

Next, consider the case when a nonzero fraction of monomers are localised in the vicinity of specific points in 
space, albeit exhibiting thermal fluctuations about these points. In this case, for many monomers (z, s) the quantity 
(5('i)(r-ci(s)))x wiU be more or less sharply peaked at some point in space and, correspondingly, (exp {ik ■ Ci{s)))x 
will not vanish identically, instead varying with k so as to reflect the spatial localisation of monomer (i, s). Then 



(exp (zk • Ci(s))) = exp (zk • bi(s)) p(.i,s)(k) 



(3.5) 



where bi(s) is the site about which monomer (z,s) is localised, and p(i^s)(k), which does not vanish identically, is 
the Fourier transform of the density of a monomer localised at the origin. In such a state, translational invariance 
is broken at the microscopic level. However, the symmetry of the state of the system at the macroscopic level is not 
settled without further information. 
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What possibilities present themselves in the situation in which a nonzero fraction of the monomers are localised? If 
the mean locations {bi(s)} of the localised monomers are distributed randomly and homogeneously over the volume 
of the system then the state is said to be macroscopically translationally invariant (MTI), the inclusion of rotational 
invariance being understood. Wc mean by this that there is no periodicity, or any other macroscopic feature capable 
of distinguishing one equilibrium state from any global translation or rotation of it. We refer to such states as 
(equilibrium) amorphous solid states. On the other hand, if the mean locations {bi(s)} of the localised monomers are 
distributed inhomogeneously over the volume of the system then the state is said to break translational invariance 
macroscopically. Examples of such states are the globular state [^7| , in which the monomers have condensed (in 
space) within a subvolume of the system, and the crystalline state, in which the mean locations of the monomers are 
arranged in a periodic lattice (and the monomers may be regarded as having condensed in wave vector space). 

How are the various possible states diagnosed by the order parameter? As we showed above, the order parameter 
is zero for all {k^, . . . ,k3} in a state that is translationally invariant at the microscopic level (i.e., a liquid). On the 
other hand, it will take at least som e nonzero values for any state in which translatio nal invariance is broken at the 
microscopic level. By using Eq. (3.5) we see that in such a state the order parameter (3.1) becomes 



1 



N 



N ^ 



y dsp(,,,)(kl)p(,,,)(k2) ...p(,^,)(k9) exp(i(kl+k2 + ...+k9) •b,;(s; 



(3.6) 



This order parameter also provides a way to distinguish between nonliquid states that are MTI and those that are 
not. In the case of an MTI state the summation of complex phase factors will totally destructively interfere unless 
the wave vectors happen to sum to zero, the random locations of the mean monomer-positions otherwise leading 
to random phase cancellations. Hence, the order parameter will only fail to vanish for values of the wave vectors 
{k^, k^, . . . , kf } that sum to zero. This property of being MTI is a fundamental characteristic of the amorphous solid 
state. In the non-MTI case total destructive interference is avoided not only if the wave vectors sum to zero but also 
under other circumstances. Hence, in this case the order parameter will also fail to vanish for certain values of the 
wave vectors {k^, k^, . . . , k^} that do not sum to zero. To establish this, consider how formula (3.6) transforms under 
a global translation by an arbitrary vector a: 



kf) •b,(s; 



[ c?sp(»,s)(k^)p(»,.)(k2)---p(»,.)(k3) exp (i(ki +k2 + • 
1=1 -^0 

^ exp (i(ki + k^ + • • • + kf) • a) 
1 1 

ds P(^,s){^^)p(^,s){'i^^) ■ ■ ■ Pi^,s){'i^') exp (^i (k^ + k^ + • • • + k^) • b, (s)J . 



(3.7) 



In situations of MTI, this transformation must leave the order parameter unchanged for all vectors a. This enforces 



the condition that for MTI situations the order parameter must vanish unless k^ 
order parameter becomes 



1 



k^H hks 



N 

TV ^ 



+ k^ H h kf = 0, and thus the 



(3.8) 



where (5pi,p2 is a d-dimensional Kronecker (5-factor, which is nonzero only if the d- vectors and have all components 
equal, in which case it has the value unity. In this state, in contrast with the crystalline state, there is no periodicity 
associated with the spatial pattern of localised monomers, and thus there will not be a collection of reciprocal lattice 
vectors for which the order parameter fails to vanish. In particular, the order parameter vanishes for 5 = 1. (One may 
equivalently regard the amorphous solid state as the special case of the crystalline state in which the unit cell of the 
crystal is the entire sample, i.e., a realisation of Schrodinger's "aperiodic solid" [p8|) The equilibrium amorphous solid 
state is characterised by the presence of random (i.e., nonperiodic) static density fluctuations, which spontaneously 
break translational symmetry at the microscopic (but not the macroscopic) level. 

In the non-MTI case, either the transformation ( |3.7| ) must leave the order parameter unchanged for a discrete lattice 
of vectors a, or it need not leave the order parameter unchanged for any value of a. When there is invariance for a 
discrete lattice of vectors (i.e., in the crystalline state), the order parameter must vanish unless k^ -f k^ + • • • + k^ = G, 
where G is any reciprocal lattice vector of the crystal (including the zero reciprocal lattice vector). When there is no 
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vector for which the invariance holds (i.e., in the globular state) the order parameter need not vanish on symmetry 
grounds for any values of the wave vectors {k"'^, . . . , }. 

To summarise, the values of the order parameter for the various values of g and the wave vectors {k^, k^, . . . , k^} 
serve to distinguish between liquid, crystalline, globular and amorphous solid states: for liquid states the order 
parameter vanishes for g = 1,2,3,...; for amorphous solid states it vanishes for all wave vectors that do not sum 
to zero (and thus vanishes for g — 1); for crystalline states it only vanishes for wave vectors that fail to sum to a 
reciprocal lattice vector (and therefore is nonzero for some {k^, k^, . . . , k^}, even when g ~ 1); and for globular states 
it need not vanish in symmetry grounds for any values of {k^, k^, . . . , k^}. 



B. A simple idealisation: generalised Einstein model 



To illustrate the general properties of the order parameter, and to motivate the specific hypothesis for the form of the 



order parameter described in Sec. IIIC and applied in Sec. VI, we examine a simple caricature of the amorphous solid 



state. We refer to this caricature as a generalised Einstein model, by analogy with the Einstein model of a crystalline 
solid adopted for the computation of the specific heat, in which it is assumed that every atom is independently 
localised by an identical harmonic potential |3^ . In the context of amorphous solidification, the caricature is obtained 
by asserting that a fraction (1 — q) of the monomers (the so-called sol fraction) are delocalised, with each monomer 
(«,s) of the remaining fraction q (the so-called gel fraction) being localised near a random mean position b,j(s), its 
location exhibiting thermal fluctuations about that mean position. We emphasise that our usage of the terms gel and 
sol in this article refers solely to the issue of whether or not a monomer is localised. Ultimately, however, we shall 
see that the gel fraction defined in this way coincides with the more common architectural definition, in the sense 
that localisation will be seen to occur only for crosslink densities for which the network spans the entire system. It 
is further asserted that the probability distribution for the fiuctuations in location of each localised monomer (i, s) 
about its mean position is gaussian and isotropic, and characterised by an inverse square localisation length Ti{s). 
Then if (i, s) is a localised monomer its Fourier-transformed density would be given by 

( exp (ik • c,(s)))^ = cxp (^k • b,(s)) exp ( - e/2n{s)), (3.9) 
so that the order parameter (prior to disorder-averaging) becomes 

(1 - 9) n Si% + ^ E / Ob.(s) • exp ( - E |klV2r.(s)), (3.10) 

a=l i=l a=l a=l 

where it is understood that the summation in the second term only includes localised monomers. 

To obtain the disorder-average of the order parameter we assume that for each monomer {i, s) the random variables 
hi{s) and Ti{s) are uncorrelated. Furthermore we assume that bi(s) is unformly distributed over the volume V and 
that the inverse square localisation length Tj(s) has the probability distribution p(t). (The notion that the state could 
be characterised by a statistical distribution of localisation lengths was introduced in pO[ | and implemented in p^.) 
In this case the order parameter becomes 



(1 



9 j-oo 9 

9)n^oi" kJ drp(r)exp(-^|klV2r). (3.11) 

a— 1 ^ a— 1 



The first term accounts for the delocalised monomers, and the second term accounts for the localised monomers. If 
q — Q then the state described by this order parameter is the liquid state. If q 7^ then it describes an amorphous solid 
state. The Kronecker 5 factor in front of the second term is a reflection of the MTI that characterises the amorphous 
solid state. This hypothesis is a reflncmcnt of the gaussian hypothesis used in a number of contexts pO|. It is useful 



to observe that the gel fraction q can be extracted from the order parameter Eq. (3.1) by taking the limit of the order 



parameter as {k^, . . . , k^} — > {0, . . . , 0} through a sequence for which X]a=i ~ ^■ 



C. Replica— order-parameter hypothesis: gel fraction and distribution of localisation lengths 



Having discussed the physical order parameter capable of diagnosing equilibrium amorphous solidification, we now 
anticipate the development of the replica approach by describing the particular form that we shall hypothesise for 
the replica order parameter, i.e., the order parameter that emerges from the application of the replica technique and 
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represents, in the replica approach, the physical order parameter discussed in the previous two subsections of the 
present section. This form is motivated by the general characterisation of amorphous sol idifica tion in terms of the gel 
fraction q and the distribution of inverse square localisation lengths p(r) given in Sec. IIIB, Below, in Sec. VI, we 
shall show that within the context of a certain model of randomly crosslinked macromolecular networks the form that 
we now hypothesise for the replica order parameter is sufficiently broad to allow us to provide an exact and physically 
appealing mean-field-level description of the transition to and properties of the equilibrium amorphous solid state of 
randomly crosslinked macromolecular networks. 



As we shall see in detail in Sec. IV, the replica representation of the physical order parameter is given by 

N 



{nY. / rf^exp(z^k".cr(.))) (3.12) 



in the replica limit, n — > 0. As we shall also see there, (• • denotes an expectation value for a pure (i.e., quenched- 
disorder-free) system of n + 1 coupled replicas of the original macromolecular system. Note the inclusion of degrees 
of freedom associated with a replica labelled by a = 0. For the sake of notational convenience we introduce hatted 
vectors {e.g., k or c), which are {n + l)(i-component vectors comprising [n + l)-fold replicated sets of d-component 
vectors (e.g., the wave vectors {k*', k^, . . . , k"} or the position vectors {c*', c^, . . . , c"}). We define the extended scalar 
product fc • c by fc • c = X]q=o ' having the special cases k'^ = k ■ k and = c ■ c. In terms of this notation, the 
order parameter becomes 

(nT. rfsexp(zfc.c.(s))) (3.13) 



By translating formula (3.11) into the replica language, through the use of Eq. ( 4. 19] ), we shall assume that the 



replica order parameter takes on values expressible in the following form: 



(1 _ q) + ^ I rf^p(^) exp(-P/2r), (3.14) 







where k = X]a=o ^ permutation-invariant d-vector built by summing the elements of the replicated vector A:, and 



^p^^'^'^^ = ria^o '^p°,q° ■ Thus, we parametrise the order parameter in terms of the gel fraction q and the distribution 
of (inverse square) localisation lengths p{t) ||41[| . This parametrisation is severely restrictive but physically plausible. 
In order for p(t) to be interpreted as a probability distribution it must be non-negative. This condition is not 
imposed a priori, but emerges from the stationarity condition. The ranging of t only over positive values reflects 
the fact that inverse square localisation lengths are positive. Moreover, delocalised monomers are accounted for by 
the term proportional to {I — q), so that p{t) must not contain a Dirac (5-function-like piece at t = 0. The factor 

iiicorporates the MTI property into the hypothesised form. It should be emphasised that the hypothesised 

form is invariant under the permutation of the replicated vectors {k'^,k^, . . . ,k"}, which is a manifestation of its 
replica-symmetric character. 



D. Symmetry properties of the order parameter hypothesis 



We now state explicitly the symmetry properties of the order parameter hypothesis (3.14) that we shall use through- 



out the remainder of this article. As we shall see in Sec. IV, the effective hamiltonian of the replica theory turns out 



to have the following symmetries: (i) independent translations or rotation of the replicas, and (ii) permutations of 
the replicas. In the liquid state the order parameter retains all these symmetries. In the amorphous solid state the 
symmetry of the order parameter is reduced. By invoking our hypothesis for the order parameter we are assuming 
that in the amorphous solid state the residual symmetries are: (i') common translations and rotations of the replicas, 
and (ii) permutations of the replicas. In other words, in the transition to the amorphous solid state the symmetry of 
independent translations and rotations of the replicas is spontaneously broken. As a consequence of the spontaneous 
breaking of certain symmetries there is a manifold of symmetry-related values of the order parameter that describe 
the solid state. 

In this article we have restricted our attention to order-parameter hypotheses that are invariant under the permu- 
tations of all (n + 1) replicas (i.e., that are replica-symmetric). This mathematical restriction is equivalent to the 
physical condition that, upon amorphous solidification, the (overwhelming fraction of the) system must exhibit one 
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member of a unique family of equilibrium states (i.e., statistical arrangements of the macromolecules) , this unique 
family of states being related by global translations and rotations. Whilst the occurrence of a unique family would not 
be an unreasonable consequence of crosslinking, especially in view of our exclusion of the anholonomic contraints that 
crosslinking introduces into the physical system, one might anticipate that crosslinking would cause the full physical 
system to exhibit many families of states (i.e., there would be states that are not related by global translations and 
rotations). Such an occurrence would be signalled by an order parameter that is no longer invariant under permu- 
tations of the replicas {i.e., for which replica symmetry is spontaneously broken). For discussions of these matters, 
see Ref. as well as Refs. [f -11 . The issue of whether or not there exists a more accurate treatment of randomly 



crosslinked macromolecular systems that would indicate the exitence of many un-symmetry-related states is a matter 
of ongoing research |^3| . 

It must, however, be emphasised that, regardless of the issue of the intactness of permutation symmetry, the primary 
physical phcnomcmon at hand in the formation of the equilibrium amorphous solid state is the spontaneous breaking of 
translation symmetry [viz., the spontaneous random localisation of macromolecules). The issue of replica-symmetry- 
breaking is not an alternative to translational-symmetry-breaking: it simply addresses whether or not a system 
with given realisation of crosslinking possesses one or many un~symmetry-related ways for the macromolecules to be 
randomly localised. To allow for the possibility that replica-symmetry-breaking accompanies translational-symmetry- 
breaking is to explore a more general class of behaviours of the system. 

E. Connection with scattering experiments 

The order parameter that we have been addressing in the present section is, in principle, accessible via neutron 
scattering experiments Q,^, at least for the case g — 2. In fact, the elastic part of the differential scattering 
cross-section (per atom) can be written as 



N dfl 



S''\q) ^ lim ( / / dsds' exp ( - iq • Cj(s; 0)) exp (iq • Ci'(s';t)) 



_^ |&incoh| dsexp(-iq-Cj(s;0)) exp(iq-Cj(s;i))^^, (3.15) 

where bcoh is the average scattering length, &incoh is the variance of the scattering length, Ci{s;t) is the position of 
monomer s on macromolecule i at time t, and (• • •) indicates a time-dependent equilibrium expectation value. The 
second part on the right hand side is the incoherent contribution, and can be extracted in some cases. By using the 
fact that the connected correlators vanish for t oo, we see that this second part reduces to 



I ^incoh I \ ^ 



N ^ 



y" ds(^exp(-iq • Ci(s)) ^ <^cxp (jq • Ci(s)) ^, (3.16) 



i.e., formula (3T) evaluated for the special case of {k^, k^, . . . , k^} — {— q, q, 0, . . . , 0}. Thus, the order parameter 
for g = 2 is proportional to the incoherent part of the elastic neutron scattering cross-section. 

Oeser et al. have measured the time persistent part of the incoherent scattering function in neutron-spin-echo 
experiments. They fit their data, which are taken in the high crosslinking limit, to a gaussian in wave vector space 
characterised by a typical length scale I, which turns out to be comparable to the radius of gyration. A potential 
critique of neutron scattering experiments results from the available time scales, of order 10 ns, which make it difficult 
to extrapolate to infinite time in order to extract the time-persistent part of the autocorrelation. This may not be 
a severe problem in the high crosslinking limit, in which one expects rather small time scales associated with small 
distances between crosslinks. However, it may become prohibitive for weakly crosslinked systems, which barely sustain 
an infinite cluster. 

Pulsed field gradient NMR (see, e.g., Ref. [^) is another experimental technique for measuring the intermediate- 
time incoherent scattering function with a spatial and temporal resolution that is complementary to neutron scattering 
experiments. Typical time scales in NMR experiments are of order milliseconds, and length scales are restricted to be 
greater than 10 nm. An example of such measurements are the detection of spatial fluctuations in swollen networks 
in Ref. M. 
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IV. REPLICA APPROACH FOR DISORDER- AVERAGED QUANTITIES 



Having prepared the way by discussing the model and the construction of a suitable probability distribution for 
the disorder, and defining an order parameter capable of diagnosing the possible states of the system, we now turn to 
the computation of disorder-averages of important physical quantities, such as the free energy, order parameter and 
certain correlators. A direct assault on this task, as it stands, seems prohibitively difficult, but it can be rendered 
tractable by the use of the replica technique [p5| , pioneered in the context of macromolecular networks by Deam 
and Edwards In this approach, we do not consider just the original degrees of freedom but, instead, a system 
comprising n+1 interacting copies {i.e., replicas) of it that will be labelled by the superscript a = 0, 1, . . . , n. In this 
new system, the quenched randomness disappears from the formulation, at the price of introducing an inter-replica 
interaction. 

A. Replica-Helmholtz free energy 

In order to compute quantities such as the disorder-averaged free energy, the order parameter and the elastic free 
energy ||4g[| in a unified way, it is particularly convenient to introduce a certain generating functional <I'„(/Lt^, {[/}), 
which we refer to as the replica-Helmholtz free energy and define as follows: 

N 



^ ^,^'^l■JO Jo 



exp(-^^d7V<i>„(A^^{^7})) ^ ^ ^ ^ ^ -. (4.1) 



n+1 



(exp(-y |] ^ ds^ ds' S^'^c^is) - c^s'))) 



The expectation value (• • is defined to be 



w 

n+1 



Yl Vc°' exp {-H, 



w ^ 

n+1 J 



/ ri 
n 2?c" exp i^H^,) 



(4.2) 



i.e., an average with respect to the [n + l)-fold replicated Wiener measure, which is proportional to exp ( — 



where 



n+l^' 



N 

I ' n ^ 

(4.3) 



Q = 4 = 1 



In analogy with magnetic systems, the properties of which one can probe by computing the free energy as a function of 
an external, position-dependent, magnetic field, it is useful in the present context to introduce the external potential 
U on which <&„ depends. The potential U , which takes as its argument an [n + l)-fold replicated d-vector x, acts 
simultaneously on replicas of all monomers, has zero average (over x), and plays the role of a source field. As we shall 
now show, knowledge of <I'„(/i^, {t^})|j/_Q allows us to compute the disorder-averaged free energy per macromolecule 
df . Subsequently, we shall see that by (functional) differentiation of <I>„(/i-^, {U}) with respect to U we shall be able 
to compute the order parameter. Furthermore, it is possible to investigate the elastic properties of the amorphous 
solid state of randomly crosslinked macromolecular networks via making a Legendre transformation with respect to 
U (see M). 
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The quantity Q^, defined by 



ds exp I ik ■ Ci{s) ) , 



(4.4) 



is the Fourier transform of the replicated monomer density. As we shall see in Sec. IVC, a certain average of Qf, 
turns out to be the replica representation of the order parameter. In order to generate Qf, it is useful to introduce 
the Fourier representation of the external potential U, viz., 



U{c) = 



V 



exp (ik ■ c) , 



(4.5) 



which leads to the replacement 



yn+l 

N 



N „i 

1 Jo 

1=1 h. 



(4.6) 



in Eq. (4.1), so that Uf. is a source for Qf,. We shall, in addition, make use of expectation values involving the weight 
that features in the replica-Helmholtz free energy, Eq. (4.1): 



• • exp 



p 

n+l,U 



xexp(-^^C/,Q,j)^^^^ 



exp 



X exp 



A2 
2 



N 



(4.7) 



ds rf,s'^5W(c?(s)-c?(.')) 

2-tr ^ /•! /•! " 

{^E ds'Y[Si^\cfis)-c?is')) 
X exp ( - ^ U^Q 



w 

n+1 

The superscript P indicates that the average is taken over a pure system; the subscript n+1 indicates that the system 
comprises n+1 coupled replicas. In fact we shall also find ourselves making use of the zero-potential expectation 
value (• • c/lt/— 0' which we introduce the special notation ( ■ • • ) 



p 



The definitions given in the present subsection give the basic mathematical objects that we shall manipulate in 
our formulation. As was anticipated above, the quenched randomness associated with the random crosslinks does 
not appear explicitly in these definitions, but instead a term proportional to /j,^ appears, which introduces a coupling 
between the replicas. In contrast with the case of conventional spin glass models the replica interaction couples all 
the replicas simultaneously, rather than pairwise. This feature is responsible for the occurrence of an order parameter 
involving a product over all the replicas, rather than pairs of replicas. For a similar feature in the context of dilute 
spin glasses, see Ref. M 



B. Application of the replica-Helmholtz free energy to the free energy / 



We now show how knowledge of 4>„(/x^, {0}) allows us to obtain the disorder-averaged free energy per macromolecule 



df, defined in Eq. (2.13). To see this, consider the quotient 

exp(-ndiV$„(^2 {0})) / exp(-ndA^$„(0,{0})). 



(4.8) 



By using Eq. (4.1) we see that this quotient is given by 
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exp(-ndjV$„(M''{0})) 
exp(-ndiV$„(0,{0})) 



N 



2N 



i,i'=l 



(4.9) 



Here, we have transformed the expectation values from ones weighted by the rephcated Wiener measure, Eq. ([4.3|), 
to ones weighted by the repHcated Edwards measure, which is proportional to exp ( — H^^i^, in which 



a=0 i=l 



n N „i „i 



(4.10) 



We indicate such expectation values by (• • the sub scri pt n + 1 indicating the presence of replicas. Next, we 

replace the exponential function in the numerator of Eq. (4.9) by its power series expansion, thus obtaining 

exp {-ndN^n (mMO})) 
exp (-ndiV$„ (0, {0})) 



2-^ Ml 

A/=0 



M! V 2N 



M N „i „i n ME 

^ ."-r. Jo Jo n ^ ' ^^+1 



^^P(W E X d.'5W(c,(.)-c,,(,s'))))^ 



(4.11) 



By expanding the M power as a multiple summation and integral over products we obtain 



exp(-ndjVj,„(^2 {o|-)^ 
exp(-ndiV$„(0,{0})) 



^ A/fl 



1 (^?V 



M\ \ IN 



N Nil 

j: ■■■ / rf^i / ds'i--- 

■ 1 ■ ; , -JO Jo 

1 1 n M g 

dsM / ds^M n n '^^'^ « (*e) - (S^)) ) 

«-p(wE / ds'<5W(c,(.)-c,(.'))))^ 



(4.12) 



We now recognise that the expectation value in the numerator factorises on the replica index to give n + 1 factors of the 
partition function Z({ie, Se; s'eleii) , Eq. ( |2.8| ). One of these factors, when taken together with the denominator 
and remaining t erms {i.e., factors other th an n partition functions), reconstructs the Deam-Edwards probability 
distribution Eq. (l2.12|) . Thus, by using Eq. (l2.14|) we obtain 



exp{-ndN<^n (m',{0})) 
exp (-ndiV$„ (0, {0})) 



ds 



ds', 



oo N N .1 .1 

E E ••• E / d.s'i... 



M 



Z({z 



Se=lJ 



By taking the logarithm of both sides, and using the expansions 



n n In z 



1 + nlnz + 0{n'^), 



ln(l + nz) ^nz + 0{n'^), 



(4.13a) 
(4.13b) 

(4.14a) 
(4.14b) 



valid for small n, we arrive at the relationship between the replica-Hclmholtz free energy $„(//^, {[/}) and the disorder- 
averaged free energy per macromolecule df of Eq. ( |2.13D : 



lim ($„(//2,{0})-<i>„(0,{0})) =/. 



(4.15) 
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C. Application of the replica- Helmholtz free energy to disorder-averaged observables 



In the previous subsection we discussed the apphcation of the rephca technique to the computation of the disorder- 
averaged free energy. We now turn to the issue of the computation of disorder-averages of equihbrium expectation 
values of physical observables as well as of disorder-averages of sums of products of equilibrium expectation values of 
physical observables. As we have seen in Sec. the latter type of quantity arises in the context of the order parameter 
for equilibrium amorphous solidification. This com putation is accomplished by using t he connection between the 
replica-Helmholtz free energy $„(/x^, {U}), Eq. (|4.l|), and expectation values of Q^, Eq. (4.4): 



—ndN<^>,,{^^^{U})^{Q,)l^^ 



+ 1,(7 



(4.16) 



More generally, by (functional) differentiation of $„(/i^,{C/}) with respect to U we can obtain average values for 
powers of the order parameter: 



ndN^n{p^,{U}) 



(-1) 



r+l 



r/n+l.. 



(4.17) 



The subscript c ind icates that the correlators are connected ones [pT| , ^ . 

Now, Eqs. ( 4.16 ) and ( 4.17 ) exhibit relationships between quantities defined in the framework of pure systems of 
replicated degrees of freedom. On the other hand, as discussed in Sec. the state of the system is appropriately 
diagnosed in terms of an order parameter built from the disorder-average of products of expectation values taken in 
the (unreplicated) physical system. The connection between the former (replica) quantities and the later (nonreplica) 
quantities is contained in the following formula, established in App. |^ which is valid in replica-symmetric states [^ : 



[(Oo)x(Oi)x • ■ • {0,)x] = lim (Oo({c^(s)}) • • • 0,({cf (.)}) 



n+l 



(4.18) 



The most important application of this connection is the case of the order parameter. For this case, as we also show 
in App. ^, Eq. ( [4.18D reduces to: 



lim — ndAr$„(^2 {[/}) 
"^0 oU? 



u=o 



hm (Qf)^ , , 



-^y^ ds(exp(zk0.c,(s)))^(exp(zki.c,(s)))^...(exp(*ks.c,(s)))^ 



(4.19) 



where 1 = {k°, k-'^, . . . , k^, 0, . . . , 0}, and {k°, k^, . . . , k^} are kept fixed as the limit n ^ is taken [|5S 
replica-limit expressions for the disorder-average of functions of the quantity (|3.1[) can also be obtained. 



Similar 



V. FIELD-THEORETIC REPRESENTATION 



The purpose of the present section is two-fold. First, in Sees. VA and VB we transform our representation of 



the physical problem from one expressed in terms of the semi-microscopic replicated macromolecular coordinates to 
one expressed in terms of generalised monomer-densities. These densities are closely relat ed t o the order parameter 
that, as we have seen, diagnoses the various physical states of the system. Then, in Sees. 



and VD 



we obtain a 

field-theoretic representation by ap plyi ng a sequence of Hubbard-Stratonovich decoupling transformations MM to the 



replica-Helmholtz free energy Eq. (5.5) This type of decoupling strategy was first used in the context of crosslinked 
macromolecular systems by Ball and Edwards ; see also Q . 

The motivation for making these transformations is that they will provide us with a suitable starting point for 
developing a mean-field description of the transition to and properties of the amorphous solid state. Furthermore, 
they provide a starting point for the analysis of fluctuations and for the investigation of the elastic properties of the 
amorphous solid state (see Ref. 0]). 
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A. Fourier representation of interactions 



At this stage it is useful to introduce the foUowing Fourier representations of the Dirac i5-functions: 



V 
p 



a— p'^',...,p^^ a— p 

where bold- face wave vectors (such as p") are d-component vectors having components taking on the usual values 
associated with periodic boundary conditions, i.e., all positive and negative integral multiples of 2tt/V^/'^ (because 
V^^"^ is the length of each side of the d-dimensional cubic container of volume V) [29|] . 

We now introduce a particularly convenient decomposition of the terms in a summation over a replicated wave 



vector p, such as that appearing in Eq. (5.1b). Consider a generic replicated vector -p = {p°,p^, . . . ,p"}. Of the 
n + 1 component d-vectors, establish the number h that arc nonzero d-vectors. Then we say that the replicated 
vector p resides in the /i-replica sector. For example, ii p = {0, 0, q^, 0, q^, 0, 0, . . . , 0} with q^ and q'' both nonzero 
d-vectors then h = 2, and we say that p resides in the 2-replica sector. The decomposition that we are introducing 
amounts to separating from the summation over p the term in the 0-replica sector (i.e., the term corresponding to 
p = = {0, 0, . . . , 0}), and also separating the terms in the 1-rephca sector {i.e., terms corresponding to those values 
of p in which exactly one d- vector is nonzero) . Thus we shall decompose summations over p into contributions from: 
(i) the 0-replica sector; (ii) the 1-replica sector; and (iii) the remainder, which we refer to as the higher-replica sector, 
and which contains the /i-replica sectors for 2 < h < n + 1. Schematically, the decomposition can be expressed in the 
following way: 

n 

E2p^2o + EE'2k + E2^- (5-2) 

where ^'y. denotes a summation over all values of the d-vector k except that the k = term is omitted {i.e., it 
comprises terms in the 1-replica sector), denotes the value of when k is in the 1-replica sector {i.e., the a'^ 
d-vector entry in k is nonzero, all other entries being zero), and denotes a summation over replicated vectors k 
residing in the higher replica sector. It will turn out to be useful for us to consider V to be larg e bu t (at least initially) 
finite, in which case it is straightforward to implement the replica-sector decomposition, Eq. ( |5.2| ). 

With a view to subsequent decoupling transformations, it is useful to use the Fourier representations of the Dirac 
(5- functions and the replica-sector decomposition in order to re-express the Dirac (5-function interactions, i.e., the 
non- Wiener measure terms that couple the replicated degrees of freedom. Thus, we see that the interaction terms in 



Eq. (4T) can be written as 

N „i 



JV „i „i n 

J2 ds ds'Y.S^'''>{cf{s)-c?{s')) 
,j'=i-^o Jo 

EE^E / dsexp{^k■cf{s}) 



a=0 k i=l 

'v 

k 



(5.3a) 
(5.3b) 
(5.3c) 

(n+l) + ^^^^|Qd'; (5.3d) 



(^+l) + TrEEbE/ rf^cxp(zk.cf(s)) 

Q=0 k 4=1 

^("+i) + ^EEUe/ rf-xp(.k.cr(.)) 



a=0 k 1=1 

-(-+!) + - 



a=0 k 

fds [\s' f[S('\cf{s)~c^{s')) 

4,4' = l-^0 -^0 a=0 
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l^v ^ 



ds exp (ifc • Ci(s)) 



(5.3e) 
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y 



V v'l ^ 

n+1 ^ 2^ hv 



TV 



a=0 k 



N ^ 



_i ^0 



ds cxp (ik • c"(s)) 



1 



N 



yn+1 A^l TV ^ 



Ar2 



TV 



T 2^ 2^ iv 2^ 



a=Q k 

+ 



1 . 2 

dsexp (ifc • Ci(s)) 



(5.3f) 



dscxp (ik • cf{s)) 



N 



]/n+l I TV ^ 

k «=1 



ds exp (iA; • Ci(s)) 



2iV 



a=0 k 



(5.3g) 



(5.3h) 



For each of the left hand sides we have performed four steps. In the first step we have used Eqs. (5.1a) and (5.1b) to 
re-expr ess t he Dirac (5-functions. In the second step we have performed the rephca-sector decomposition, according 
to Eq. (5.2). In the third step, we have recognised that the summands in the summations over wave vectors of the 
second step are even functions of the relevant wave vector. Furthermore, none of the summations includes a zero wave 
vector. Thus, in each case the summation can be restricted to half of the relevant wave vector space, provided a factor 
of two is included to compensate. To represent this, we have introduce the notation to denote but with k 

restricted to the half space via the additional condition k • n > for a suitable unit c?-vector n, and to denote 
y^j^ but with k restricted to the half space via the additional condition k ■ h > for a suitable unit (n + l)(i- vector 
n. The virtue of this procedure is that in our subsequent development it will enable us to avoid the introducti on o f 
kinematically non-independent fields. In the fourth and final step we have used the definition of given in Eq. ( [4.4| ), 
which in the context of the replica-sector decomposition becomes 



Qk = 



1 ^ /•! . . 

— / ds exp f zk • cf (s) j , (1-replica sector), 

— J ds exp yik ■ Ci{s)j , (higher-replica sector). 



(5.4a) 



(5.4b) 



B. Applications of the Fourier representation 

By applying the results of Sec. |V A| to Eq. ( [4.l[ ) we obtain 
exp(-nd7V$„(/x2,{;7})) = 



a=0 k 



exp 



A§Af2^t|„n|2\\W 



V 

k 



ri«i)). 



where we have introduced the effective excluded-volume parameter 



(5.5) 
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i.e., the bare excluded- volume parameter renormalised to a smaller value by a correction term proportional to the 
crosslink density parameter /i^. The prefactor Bn is an unimportant constant, which arises from terms in the 0- replica 
sector, and is given by 



6, 



We see from Eqs. ( |5.5|) and (5^) that in the 1-replica sector there is a competition between the excluded- volume 
interaction and the effect of the crosslinking, represented by /j,^. If A^ is positive then configurations having 
nonzero are disfavoured (and therefore MTI is favoured), whilst if A^ is negative then configurations having 
nonzero are favoured (and therefore there is a tendency to violate MTI). As we are primarily concerned with 
investigating the liquid and the amorphous solid states, both of which are MTI, we focus our attention on the regime 
A^j > 0. On the other hand, the (higher-replica sector) component of the term due to crosslinking increases the 



statistical weight of configurations in which Qf. ^ (for k in the higher-replica sector). As we shall see in Sec. VIB, 
the coefficient of this term, /i^, is the control parameter governing the transition from the liquid state (for small /i'^) 
to the amorphous solid state (for large /x^), as characterised by the values of the order parameters shown in Table |. 



C. Hubbard-Stratonovich decoupling scheme: replica-Helmholtz free energy at zero external potential 

We now obtain a field-theoretic representation by a pply ing a sequence of Hubbard-Stratonovich decoupling trans- 



formations to the replica-Helmholtz free energy, Eq. (5.5). We focus on the case J7 = in the present subsection, 



and in the following subsection present results that allow for nonzero U. The motivation for these transformations 
is that via them all interactions between different macromolecules are eliminated at the expense of introducing a 
certain additional stochastic field to which the monomers are coupled. This strategy has the following virtues. First, 
the task of summing over the configurations of the system of N replicated macromolecules is reduced to the task of 
summing over the configurations of a single replicated macromolecule, albeit one that is coupled to the stochastic 
field. The monomers that constitute this replicated macromolecule remain coupled to each other via the Wiener 



measure, Eq. (2.2), and b y the stochastic field to which they are coupled. To the extent that this summation can 



be performed, Eq. (5.12a) gives the replica-Helmholtz free energy in terms of functional integrals over the sto chast ic 



fields. Second, the stochastic field itself has a natural physical interpretation: as we shall see explicitly in Sec. VP it 
is related in a direct manner to the order parameter. 

The appropriate Hubbard-Stratonovich decoupling transformations are predicated on the multiple use of the fol- 
lowing pair of integrals: 

exp ^ — a|u;|^^ = (o/tt) J dRe z dim z exp ^ — a|z|^^ exp ^2m Re zw*^ , (5.8a) 
exp (^ + a\'w\^^ = (a/n) J dRe z dim z exp ^ — a|z|^^ exp ^2a Re zuj*^ , (5.8b) 

where w is an arbitrary complex number, a is a real and positive (but otherwise arbitrary) number, and the inte gral s 
are taken over the entire complex z plane. We transform each exponential term inside an expectation value in Eq. (5.5) 
(with J7 = 0) by using these integrals, those having the coefficient A^ or Aq with Eq. ( ^.8a ) and those having the 
coefficient /i^ with Eq. ( 5.8h| ). 



We define the measures V}-^, V and via 



" 4. \ 2 at2 

Vln = n n -^^^^^^ dlmf]£, (5.9a) 

Q = k 

V^'n = t^dReflj^ dlm%, (5.9b) 



fe 



Vloj = TT -2 dRewk dimwk, (5.9c) 

k 

where Ylk denotes the product over all d- vectors k in the half-space given by the condition k • n > for a suitable unit 

d-vector n, and hI denotes the product over all (n -I- l)d-vectors k in the half-space given by the condition k ■ n > 
for a suitable unit {n + l)d- vector h. These definitions have convenient normalisation properties [^5[: 
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(5.10a) 



/ V^'ncxp ( - tJ-'^NV-''J2i,\^k\^) = 1' (5.10b) 
J Pjwexp Ag7VV~i^|^|wk|^) =1. (5.10c) 



By using this sequence of Hubbard-Stratonovich transformations Eq. (5.5) becomes 
exp(-ndAr$„(^2 {0})) = 

^ exp (2*A2 TV ^ ^^Re QZ + 2^,^NV--Y.k^'' ^0 /„+i 



J Pjc^exp ( - ~XIN^V-'Y!J^^0 { {'^^XIN^V-'Y.IrccuI Qi) 



(5.11) 



By examining the expectation values (• • and (• • in Eq. ( ^.ll| ) we see that indeed the sequence of Hubbard- 
Stratonovich transformations has led to the decoupling of the N (replicated) macromolecules from each other. More- 
over, as these expectation values are products of identical factors, one for each replicated macromolecule, we may 
replace them by the N^^ power of an expectation value involving a single replicated macromolecule. Thus, we see 
that Eq. is given by the quotient of partition functions of stochastic fields. 



Vln f v''nexp{-ndNJ'n{{^k:^k})) 

exp ( - ndN^n{f^\{0})) = Bn ^- , (5.12a) 

/ Pjwexp (-ndiVJ^'i™({wk})) 

governed by the effective hamiltonians JF„({r2jJ, fi^,}) and JF^™({ti'k}), where 

" t 1 

Q = 

-ln<Jexp(2iA2w-i^^^Rel];j* / ds exp (ik • c"(s)) ) 

xexp(^2^V""^.Rer2|y ds exp (iA: • c(s)) j ^ ^, (5.12b) 

n(ij^'*™({tJk}) =A§iVl/-i^J^|tJkf -ln{Jexp {2iXlNV^^Y^YieLulj^ dsexp (ik • c(s)) )^^. (5.12c) 



We can compute ^„({rig, ^j.}) and .F'^°"({a;k}) perturbatively |Q, order by order, in powers of {Og, fi^} and {wk}- 
This perturbative construction of the Landau-Ginzburg- Wilson effective hamiltonian is equivalent to that arising in 
many other contexts in statistical physics. 



D. Hubbard-Stratonovich decoupling scheme: repIica-Helmholtz free energy at nonzero external potential 

Wc now turn to the general case of the replica-Hclmholtz free energy, in which U need not vanish. As with the 
f/ = case, it will be convenient to perform a Hubbard-Stratonovich decoupling transformation on $„(/i^, {U}). This 
transformatio n is p redica ted o n the multiple use of the following pair of integrals, which are generalisations of those 



given in Eqs. (5.8a) and (5.8b) 
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exp 



exp 



a^wl^ - jw - jw*)^ ^ (a/7r) exp (ajj) 
X J dRe z dim z exp — a\z\^ + 2iaRe zw* — ia {zj + z*j)^, 
(a(|w|^ -jw-jw*)'j = (a/7r)exp ( - ajj) 

X J dRe z dim z exp ^ — o 1 2; | ^ + 2a Re zw* — a (^zj + z*fj^, 



(5.13a) 



(5.13b) 



where w, j and j are arbitrary complex numbers, a is a real and positive (but otherwise arbitrary) number, and the 
integrals are taken over the entire complex z plane. Following a strategy analogous to that used above for the U = 
case, we find 



ndN^nip"^, {[/}) = - lnS„ + In J Pjwexp {^ndNT'^"'{{LUk})) 



In / Vinv'nexp ( - ndNTn{{^l, %}) + * E E'^^k^^k - Yl^k % 



a=0 k 



(5.14) 



where the effective hamiltonian Tn is given in Eq. (5.12b). Thus we have transformed our description into an effective 
one for the fields {^1^' ^fci- should be noted that the coefficients of the powers of {U^, C/^} in the functional Taylor 
series for $„(/i^, {[/}) are simply related to the connected correlators of {i^^' ^fcl- 



V 



yn 



nd7V$„(M^ {f/}) = ndiV$„(/i2, {0}) - ^ E'^k C/^k ^ ^ 



E 



r(-i) 



r,s— 
(r,s)/(0,0) 



EE---EE E---E^k----f/k:f/^.---f^. 

Ql=0 ki Q,. = k^ fcj fc^ 



r! s 



n+l,c 



(5.15a) 



where the expectation value (• • •)^_|_]^ is defined via 



ri+l 



y" VlnV^'n • • • exp - ndA^J^„({17£, 17^}) 
pl^llP^f} exp ( - ndNTn{{^k,^k}) 



(5.15b) 



The virtue of the present development is that it allows us to construct quantities of physical interest, which 
typically involve the densities {Q^, Qj^} in terms of more readily computable quantities, which involve the stochastic 
fields {fl^,V,f,}. Indeed, by using Eqs. (4.17) and (5.15a) we see that 



(5.16a) 
(5.16b) 



Thus, by using Eq. ( [4.19]) we see that we can relate the order parameter to the expectation value of the stochastic 



field flf^: 



/ ds(exp(*k0.c,(.s)))^(exp(*ki.c,(s)))^...(exp(*kf .c,(s))), 
i=i "^0 



= lim (ni)^^, 



(5.17) 



where I = {k°, k^, . . . , k^, 0, . . . , 0}. In the following section we make explicit use of this development in order to 
compute the order parameter. 
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VI. SADDLE-POINT APPROXIMATION IN THE CRITICAL REGIME 



In the preceding sections we have developed an exact, formal, field-theoretic representation of the statistical me- 
chanics of randomly crosslinked macromolecular networks. In the present section we shall explore the properties of 
such systems, focusing our attention on the regime of crosslink densities near to the equilibrium phase transition from 
the liquid state to the amorphous solid state that sufficient crosslinking causes. We shall do this by analysing the 
field-theoretic representation at the level of mean -field theory, considering in detail expressions for the free energy and 
the order parameter. Following this, in Sec. VII, we shall consider the implications of a certain, physically important 
class of fluctuations. In order to streamline the presentation, a considerable amount of technical detail has been 
relegated to appendices. The main results of this section have been briefly reported in Ref. | |T^ . 



A. Approximation strategy 



The mean-field level of approximation follows from computing the functional integral in the numerator of Eq. (5.12a) 
by using the saddle-point method. This amounts to replacing the functional integral by the value of its integrand that 
is stationary with respect to variations of {^^, ^fc}? so that, omitting unimportant constants, we obtain the following 
approximations for <i>„(^^, {0}), for {^^)n+i, for (^fe)n+i and for {Q)n^i: 



$„(M',{0})=^„({f^£,f^fe}), 

*('9k)n+l = (^k)n+l — J 

where fl^ and f2^ make JF„({S7^, f2^}) stationary, i.e., satisfy the stationarity conditions 



sn*. 



0, 



= 0. 



By using Eq. ( 5.12bD the stationarity conditions become 

' / (itexp(iq-c"(i)) exp(2iA^iVy"^^^^Rei7£* / ds exp (ik • c"(s)) 

'Jo ^^g Jo 



X exp (^2/xV""^ .Re f^i / ds exp {ik ■ c(s)) 



w 



n+l 



/ " t _ 

exp hiXlNV-^ y^Yl R-cf^r / exp(ik-c"(s)) 

/ — ^ - \ \ 

X exp (^2/i2y-«^_Rer2t / ds exp {ik ■ c{s))jj 



/ dtexp(iq-c{t)) exp (2iXlNV-^y2yi ^^^k* / ds exp (ik ■ c" {s))) 



exp (2/^V-"^.Rerji / ds exp {ik ■ c(s)) 



w 



n+l 



/ " t _ 

exp 2i~XlNV-^ ^® ^k* / ds exp {ik ■ c" 

/ 1 _ /-i 

X exp (^2;u2y-"^.Re 171 / ds exp {ik ■ c{s) 



w 

n+l 



(6.1a) 
(6.1b) 
(6.1c) 

(6.2a) 
(6.2b) 



(6.3a) 



(6.3b) 
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B. Instability of the liquid state 



In the context of the mean-field approximation, the hquid state corresponds to f2g — flj^ — 0, which can readily be 
checked to solve Eqs. ( 6.3a ) and ( 6.3b ). To address the stabihty of this state we examine ^„({rig, f^j.}), Eq. ( 5.12b ), 
and expand perturbatively about fi^ = 17^, = to second order in fi^ and Sl^,. This gives 



nd^4{nZ,nf^})^\lNV-'J2J2A^ + '^nNV-'goi\k\^)) 



nr. 



(6.4) 



The correlators necessary to calculate the terms in this expansion are computed in App. and the function 5o(|kp) 
resulting from the subsequent arclength integrations is defined in App. ^ and has the value 



90 [ 



(l-iP) 



1 (i^^) 



J l-P/6, if fc2 < 1; 
\4/fc2 iffc2>l. 



(6.5) 



As we anticipated at the end of Sec. |VB| , the stability of the 1-replica sector is controlled by the coefficient of the 
[ri^l term in this expansion. This coefficient, together with the positive-definiteness of godkp), show that provided 
the crosslink-renormalised excluded- volume parameter A^j, Eq. ( |5.6| ), is positive, the 1-replica sector is locally stable. 
Thus, the saddle-point value fig is zero. The positive-definiteness of requires that 



A^ > 



(6.6) 



i.e., that the repulsive character of the physical excluded- volume parameter A^ is sufficiently strong to enable the 
system to withstand the effective tendency towards collapse afforded by the crosslinking. Thus, we see that even 
at the level of mean-field theory it is only as a consequence of the presence of the excluded- volume interaction that 
the system can, at the same time, be stable with respect to collapse to the (inhomogeneous) globular state and yet 
unstable with respect to the formation of the (macroscopically homogeneous) amorphous solid state. 



The stability of the higher-replica sector is controlled by the coefficient of the | il^ | 
i.e., by 



term in the expansion, Eq. (p.4|), 



l-^l'v-goim 



(6.7) 



(considering, as we do, > 0). The two contributions to this coefficient enter with competing signs, owing to 
the attractive nature of the effective term arising from crosslinking, and thus provide the opportunity for the loss 
of positivity of this coefficient. Indeed, the coefficient indicates that the liquid state will be stable for /i^ < 1 and 
unstable for /i^ > 1, i.e., stable only for sufficiently small crosslink density, the factor of in Eq. (3.7) being 
eliminated by first taking the limit n 0, and subsequently taking the thermodynamic limit (V oo, N ^ oo, N/V 
fixed, /i^ fixed) ||5^. The least stable modes correspond to long wavelengths, 0, for which god^P) ~^ 1 from 

below Hi. 

The linear stability analysis of the present subsection indicates that the liquid state, as characterised by the order 
parameter discussed in Sec. [II, is stable when the mean number of crosslinks per macromolecule [M]/N is smaller 
than a certain critical value M^/N , i.e., those mean crosslink densities corresponding to fi^ < 1. However, for larger 
crosslink densities, {[M]/N) > [M^/N), i.e., fi^ > 1, the liquid state is unstable j5^, being replaced by an alternative 
state which, as we shall see in the following two subsections, is an amorphous solid state, characterised by fij, ^ but 
= 0. In fact, the state that replaces the liquid state will turn out to have the property of macroscopic translational 



invariance (see Sec. ill A), so that even though it has ^ this is compatible with and does not disturb the fact 
that the 1-replica sector remains stable and that fig remains zero. 



C. Free energy 



We now set about exploring the nature of the amorphous solid state with respect to which the liquid state is 
unstable for fi^ > 1. Initially, we do this by following the strategy outlined in Sec. VIA of making stationary with 
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respect to the fields D,^ and flj^. However, we are unable to parametrise the entir e spa ce of pos sible fields. Instead 



we consider the class of fields for which physical motivation was presented in Sec. [II B [see Eq. ( 3.14 )], 

— 0, (1-rephca sector), (6.8a) 

^k^l^'^o J rfTp(r) exp ^— P/2Tj , (higher-replica sector), (6.8b) 

evaluate J-'n for such fields, and make the resulting quantity stationary with respect to the variational quantities, the 
gel fraction q (a number) and the distribution of (inverse square) localisation lengths p(t) (a normalised function). 
This amounts to making a variational mean-field approximation. However, as we shall see in the following subsection, 
the hypothesis we make for the saddle point will actually turn out to contain an exact saddle point of !Fn- 



By inspecting Eq. (5.12b) and employing Eq. (6.8a) we see that there are two contributions to the free energy: a 
term quadratic in fl^, , and a term that can be identified as the logarithm of the partition function of a single re plica ted 



macromolecule coupled to Qj^ . The explicit details of the evaluation of these terms when fi^, is given by Eq. ( p.8h ) is 
presented in App. Then the variational mean-field approximation to the free energy is given by 

/ = hm $„(/i2,{0}) « lim min .F„({1)C, %}) « min /™'{g,p}, (6.9) 

ri— >0 n^0n°,Oj; 9.p(t) 

where we have omitted constants independent of the variational parameters q and p(t), and the variational free energy 
f™{q,p} is given by 

r'{qM = -i(cxp ( - ^,\) - (1 - - In {v'/''/27:c 

,2r„r /•! 



I f°° f°° / 

+ -M'<z'y^ dTip{Ti) dT2p(r2)ln(^(Tfi+r2-i) 
2 7^1 Jo 

/•oo 

xj drlp(Tl)•••dT,p(r,)ln(w('')DetM7^J,;^,j . (6.10) 

Here, Det'-'^' denotes the determinant of an r x r matrix, Ti-^ppi is an (r x r)-matrix-valued function of the r arclength 
coordinates {si^j^^i and the r inverse square localisation lengths {t^Yv=i^ ^-nd W'^'^^ is a single such function, Ti-^pji 
and W*-''^ being respectively defined in Eqs. (Fl) and ( F2b ) of App. ^ 



As anticipated in Sec. [IF, in addition to intensive terms we find a nonintensive term, proportional to InV^, owing 
to the omission of the disorder-average of the Gibbs symmetry factor. The presence of the In V factor signals the fact 
that the configuration integral produces additional powers of V . These powers of V can only be due to degrees of 
freedom that are allowed to vary over the entire sample, i.e., to the fraction of macromolecules that are delocalised. 
Thus this term depends on the gel fraction g, but it cannot (and does not) depend on p(t), which only describes the 
localised degrees of freedom. In the following subsection we shall analyse the self-consistency condition for the order 
parameter directly and, although no quantity proportional to In V will appear, we will re-obtain the exact same results 
as in the present subsection. This approach will be see n to have the additional substantial virtue of demonstrating 



that the hypothesised form of the order parameter, Eq. (3.14), used in the present section as a variational hypothesis 



in fact provides an exact saddle point of the free energy, not merely a variational approximation. 

As a first step towards minimising /™'' we regard the term proportional to In V as dominant, and make it stationary 
with respect to the gel fraction q. This leads to the condition pO| , ^ , p^ 

exp(-^\) = 1-g, (6.11) 

For all values of /i^ this equation has the root g = 0, corresponding to the liquid state. However, for /^^ > 1 an 
additional root appears, emerging continuously from g = at /i^ = 1, and describing the equilibrium amorphous solid 
state. In Fig. |l| we show the dependence of the gel fraction on /i^. For /i^ 3> 1, i.e., the highly crosslinked regime, q 
approaches unity asymptotically as q ^ 1 — cxp ( — /i^). In the critical regime, < /i^ — 1 ^ 1, it is convenient to 
exchange /i^ for the new control parameter e, defined via 

/i^ EE 1 -t- e/3, (6.12) 



with < e <C 1. We may then solve perturbatively for g, obtaining 
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g = 2e/3 + 0(e2). 



(6.13) 



Having determined the condition satisfied by q we now turn our attention to the dependence of f™' on the distri- 
bution p{t) of inverse square locahsation lengths. As we are primarily interested in crosslink densities in the vicinity 
of the vulcanisation transition (i.e., < e <C 1), w e use the result that, to order e, we have q = 2e/3. This allows 
us to retain in the summation over r in Eq. ( |6.10 ) only the terms r = 2,3 (the r — 1 term vanishing identically by 



the construction of W*- ''■'). Next, we assume that the inverse square localisation lengths having appreciable statistical 
weight in p(t) are also of order e, i.e., small compared to unity, in units such that IL/d = 1, so that localisation 
is on length scales much larger than the size of a free macromolecule. (We shall confirm the consistency of this 
assumption a posteriori.) Thus, we may use the result from App. ^ to expand In (W^'"' Det'-'^-' T^^pjr) in Eq- (|6.10| ) for 
small {Tp}'p^i, retaining terms to order Tp. Then we integrate over the arclength variables {si, . . . , Sr} by using the 
results of App. ^. Omitting terms that are independent of p(r) we find that, to 0{e^), 

l/2e\^r fTi+T2\) 1 f 2e\^ ( f ti + T2 + 1 /2e\^f T1T2 



where the curly brace carrying the subscript r indicates averaging over the localisation lengths, i.e., {^'(ti, T2, . . .)}^ = 
dTip{Ti)dT2p{T2) ■ ■ ■ ^'(ri, r2, . . .), as many lengths ri, T2, . . . as feature as arguments of the arbitrary function . 
For future reference we note that if we suppose that the distribution of localisation lengths is sharp, i.e., has no 
fluctuations, so that p{t) = 5{t — ^~^), where £, is the sharp value of the localisation lengths, then the expression for 
yvar gjniplifies, becoming 

1 /2e^^ ,^0, 1 f2eV 1 



correct to 0{e^). In this case, demanding that /™' be stationary with respect to yields, to 0{e), 

1 _ / 0, if ^2 < J 



C^-\2e/3, ifp'>l. (^-1^) 

We now return to the general situation, in which the distribution of localisation lengths is not constrained to be 
sharp. Rather than demand that f™ be explicitly stationary with respect to p{t) itself, it is convenient to exchange 
its dependence on p{t) for dependence on the Laplace transform p{f) given by 

POO 

p{t) = / c?Tp(r) exp(— fr). (6-17) 
"'0 

The details of this exchange are deferred to App. what results is the following expression for /™', correct to 0{e'): 



1 /2e^ ^ 
12 



df{dp/dfY. (6.18) 



This expression has the virtue of being a local functional of p(f ), so that the consequent stationarity condition will be 
a differential equation for p{f ). Moreover, the (global) constraint that p{t) be normalised to unity, dTp{T) = 1, is 
exchanged for the (local) boundary condition p(0) = 1. 

We now demand that f™ be stationary with respect to p{t), i.e., that 5 / 5p{t) = 0. The details of computing 
the functional derivative of /™'' with respect to p{f) are deferred to App. what results is the following stationarity 
condition, correct to 0{e^): 

Q_^/™' /^2e\^ 1 ^^^^^\ , f2e\^ 1 ^,^,2n (2e\^l<Pp 

5p{f) 

or, equivalently, 



(l)i('-«^))-(l)i(-«^)v(|)iS, 
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f 



^ = 6p(f)(l-p(f)), (6.20) 



correct to 0(6). Normalisation of p(t) leads to the boundary condition p{Q) = 1. As p(t) does not contain a (5-function 
contribution at r = (see [^), p obeys the additional boundary condition p{oo) — 0. 

Before solving the stationarity condition we note that p(r) depends parametrically on the crosslink density, so it 
would be more accurate to denote it by p(r; e). We now introduce the scaling function it{9) in terms of which p(r; e) 
is given by p{t; e) = (2/e)7r(2r/e). In other words we transform the dependent and independent variables as follows: 

epiT;e)/2^7r{e), (6.21a) 
T^e9/2.. (6.21b) 

In this way, up to an elementary factor, the dependence of p(t; e) on t and e is combined into a dependence on a 
single scaling variable 9; see [p3[. Then the Laplace transform of the scaling function tt{9) is defined via 



if (9)= d9Tr{9)e^p{^99), (6.22) 
Jo 

so that 

p{f) = tt{9) (6.23a) 

ef/2 = 9. (6.23b) 

In terms of 7r(0), and neglecting 0{e) contributions, the stationarity condition then becomes 

9^ = 2m[l-n{9)), (6.24) 

subject to the boundary conditions tt{0) — 1 and 7r(oo) — 0. 

We have been unable to solve this nonlinear ordinary differential equation for 7r(0) analytically. One might consider 
solving this differential equation numerically, and then inverting the solution numerically to obtain t:{9) and hence 
p{t). Whilst this is possible in principle, the numerical inversion of Laplace transforms is notoriously unstable. 
Instead, we have found it preferable to take the inverse-transform of the differential equation analytically, and thus 
we obtain the following nonlinear integro-differential equation and constraint for n{9): 

= (1 - ^) A0) I' d9'nie')ni9 - 9'), (6.25a) 

d9Tr{9) = l, (6.25b) 

the constraint resulting from normalisation. 

We shall obtain 7t{9) [and hence p(t)] in Sec. VIE, and discuss the consequences of the physical values of q and 
p{t). Before doing so, we shall adopt a different point of view, in which we focus not on the variational extremisation 
of the free energy but instead on the self-consistent equation for the order parameter. 



D. Self-consistency condition for the order parameter 



In Sec. VI C we enforced the stationarity of the effective hamiltonian only with respect to the parameters q and p{t) 
of our order parameter hypothesis, and not with respect to arbitrary variations. As a consequence, we are not yet in a 
position to address whether or not the resulting order parameter is a true saddle point of the effective hamiltonian. In 
the present subsection we establish that the solution that we have found is indeed a true saddle point of the effective 
hamiltonian by directly analysing the stationarity conditions ( 6.3a ) and (6.3b) themselves. We emphasise that this 
approach allows us to circumvent the difficulties discussed in Sec. II F that arise in the computation of the contribution 
to the free energy associated with changes in the indistinguishability factors introduced by the crosslinks. The results 
of this subsection have been briefly reported in Ref. [ p^. 

We i nsert the hypothesis given in Eqs. (6.8a) and (6.8b) into the stationarity conditions (6.3a) and ( 6.3b| ) derived 
in Sec. VIA to obtain 
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(l-g)J-Q + 9%o / drpir) exp{-f/2T) 
Jo 

dt exp (^iq ■ c(t)j 



exp (^^V""(7^.(5k,o / dTp(T)exp(-P/2r) / dse 



1 - \ \ w 

ik-c{s) 



n+1 



w 

n+l 



(6.26) 



Here, in both the numerator and the denominator we have relaxed the constraints on the summations having coefficient 
fj? by: (i) doubhng the range of the summations to include the entire higher-replica sector by making use of the 
property of the hypothesis 51^ = 0*-; (ii) including the 1-replica sector terms (which vanish by the MTI of the order 

parameter hyp othesis); and (iii) inserting identical factors associated with the 0-replica sector. It should be noted 
that Eq. ( |6.26 ) also follows from the direct application of the Weiss molecular-field method. 

As shown at the end of App. ||, evaluation of the numerator and denominator of the right-hand side yields 



i^- + qSq^o / dTp{T)e 



Jo ^- Jq 



X / dTi---dTrp{Ti)---p{Tr)S{T-'T^''^), 
JO 

Sr+l,r+l - TTTm'^^P^^ Sp,r+1 " X! ^r+'^'P'^PP' ^P' -T+l^ 



(6.27a) 
(6.27b) 



p=i 



p.p'=i 



where the limit n — > has be en t ak en ev er ywhe re exc ept i n the dimension of q, and where Sppi , U^"^^ , W^'"' and C^-^^ are, 
respectively, defined in Eqs. (BS). (F2a), (F2b) and (F2c), and depend on {si, . . . , s^+i} and {ri, . . . , tv+i}. It should 
be emphasised that Eq. (3.27a) is not solved by any sharp distribution of localisation lengths p(r) = 5(r — ^~^). Thus, 



a variational hypothesis involving a sharp distribution gives at best a variational approximation, whereas a variational 
hypothesis involving a non-sharp distribution has the potential to yield an exact saddle point, and we shall find such 
an exact saddle point below, at least in the vicinity of the vulcan isation transition. 

We now extract information about q and p{t) from Eq. ( 3.27a ). First, we take the limit q^ — > 0, via a sequence for 
which q = 0. In this limit, the left hand side becomes q, and on the right side each integral gives a factor of unity, 
yielding the self-consistency condition for the gel fraction q: 



(6.28) 



i.e., precisely the self-consistency condition for q found from the free energy approach is Sec. VI C and discussed there. 

Having decoupled the is sue of the gel fraction q from the issue of the distribution p{t) we now return to p{t) 
itself. By considering Eq. (6.27b) for a fixed nonzero value of q^, and using Lerch's uniqueness theorem for Laplace 
transforms we find that indeed the hypothesis solves the self-consistency condition for provided that 

the distribution p{t) satisfies the condition 



OO 9r r /•! 



dsi ■ ■ ■ dSr+l / dri ■ ■ ■ dTrP^Ti) ■ ■ -piTr) 6 [t — Tr) 

Jo 



(6.29) 



This equation for p{t) is, for all values of fi^, identically satisfied if g = 0. 

We have not, thus far, made any approximations beyond that of mean field theory. In order to render Eq. ( 6.29| ) 
tractable, we now re strict our attention to the vicinity of the transition regime in the solid state, i.e., to values of 
e, as defined in Eq. ( |6.12 ), satisfying < e <C 1. This restriction allows us to assume that q is small, and that only 
localisation lengths much larger than the free-macromolecule radius of gyration have an appreciable probability, i.e., 
p{t) only gives appreciable weight for < r <§C 1. Thus, we need retain in Eq. (6.29) only terms for which r is 1 or 
2, and may expand T*^^^ to 0{t^) and T^^^^ to 0{t^). As discussed in App. ^, in terms of the scaling function tt{9) 
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introduced in Eq. (6.21a), we recover Eq. ( |6.25a| ) subject to the normalisation condition Eq. (6.25b), i.e., precisely 
the stationarity condition for p(r) found from the free energy approach. 

Thus, the condition that the order parameter be self-consistent turns out to be identical to the condition that the 
effective hamiltonian be stationary with respect to variations within the subsp ace sp anned by the hypothesised form 
of the order parameter. The form for the order parameter hypothesised i n Eg. (3.14) is not merely a variational form 
but in fact gives an exact saddle point of the effective hamilton ian, E q. ( 5.12b ). 

We have obtain ed th e equa tion for the gel fraction g, Eq. ( |6.11 ), and the equations for the scaled distribution 
7r(0), Eqs. (3.25a) an d (335b| ), from two different points of view. In the previous subsection we have discussed the 
consequences of Eq. ( |6.11 ) for q. In the following subsection we shall discuss the solution of Eqs. (6.25a) and ( 6.25b| ) 
for 7r(0), and elaborate on the physical consequences of our results for q and p{t). 



E. Characteristics of the amorphous solid state 



For the sake of completeness we first restate the results concerning the gel fraction q that were found in Sec. |VI C 
from the self-consistency condition on q, Eq. (6.11). For all values of we find the solution q = 0, corresponding 



to the liquid state. For > 1 an additional solution appears, emerging continuously from g = at /i = 1, and 
describing the equilibrium amorphous solid state, as shown in Fig. |^. For /i^ 3> 1, i.e., the highly crosslinked regime, 
q approaches unity asymptotically as g ~ 1 — exp ( — /i^). In terms of the deviation of the crosslink density from 
criticality, i.e., e defined in Eq. ( |6.12 ), the critical regime is < e ^ 1. In this regime we may solve perturbatively 
for q, obtaining q = 2e/3 + 0{e^), Eq. ( |6T3| ). 

It should be noted that the stationarity condition on q, Eq. (6.11), is precisely the condition obtained by Erdos 
and Renyi in the context of random graph theory |]60| , which can also be interpreted as a mean-field treatment of 
percolation. In particular, Erdos and Renyi showed that that for a random graph of TV points and ij^N/2 edges the 
probability for the fraction of points in the largest component to differ from the solution q of Eq. ( |6.1l| ) vani shes in 
the N ^ oo limit. A related approach to the theory of macromolecular networks |Q has also led to Eq. ( |6.11 ). This 
is physically quite reasonable: one would anticipate that the transition from liquid to solid would occur when the 
density of crosslinks is sufficient to create a macroscopically extended network of crosslinked macromolecules. 

We now address th e distr ibution of localisation lengths via th e scaling form 7r(£). We have solved both the intce ro- 
differential equation ( 3.25a ) and the differential equation ( 6.24 ) numerically |^,^, and the solution of Eq. ( 3.25a ) is 
shown in Fig. ^. As we see in this figure, the scaling function Tr{9) has a single maximum near = 1, away from which 
it decays rapidly. In fact, states for which Tr{6) takes negative va lues ar e not ruled out by the hypothesis Eq. ( [3.14[ ), 
but are not found as solutions of the stationarity condition Eqs. (6.25a) and ( 6.25bD . 

We are able to obtain asymptot ic pro perties of tt{9) analytically. The asymptotic form n{6) ^ aO^^ exp ( — 2/6') 
(for <C 1) is obtained from Eq. (6.25a) by neglecting the second term on the right hand side. Notice the essential 
singularity at the origin: tt{6) vanishes very rapidly indeed as 6* ^ 0. The coefhcient a ~ 4.554 cannot be obtained 
from local asymptotic analysis. Instead we have obtained it separately by the numerical solution of Eq. ( |6.24 ), as 
discussed in footnote The asymptotic form 7r(0) ^ 3(b6 — 3/5) exp ( — 66*) (for 6* ^ 1) is obtained by computing 

the inverse Laplace transform of the approximate analytical solution of Eq. ( |6.24 ) near the point 9 = —5, at which 
7t{6) diverges. The coefficient b ^ 1 .678 was obtained separately by determining the (negative) value of 9 at which 
the numerical solution of Eq. ( 6.24 ) diverges. Notice the exponential decay of Tr{9) for large 9: Tr{9) goes to zero 
quickly as 9 —^ oo. For the sake of comparison with the numerical results, the small- and large-0 asymptotic forms 
for 7r{9) are also shown in Fig. ||. A distribution of localisation lengths also features in Panyukov's approach to the 
well-crosslinked regime; see ]54| . 

In addition to providing the distribution of localisation lengths, knowledge of tt{9) allows us to construct the order 



parameter, fi^,. By using Eqs. (6.21a) and (6.21b) in Eq. ( |3.14 ) we obtain 



% = (l-2e 



/3)<5^a-f (2e/3)<5^_oC.(y2^), 
(fc)= / d9T:{9)e-^"'^'. 



(6.30a) 
(6.30b) 



Although we do not have an exact analytical expres sion fo r 7r(^), we can compute Lo[k) numerically. We do this by 
inserting the the numerical values of -k{9) into Eq. ( |6.30b ), and show the result for uj{k) in Fig. ^. We are able to 
obtain analytical asym ptotic expressions for u}(k). For q <C 1 the result simply follows from expanding the exponential 
function in Eq. (3.30b) in powers of g, thus obtaining 
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Uj{k) 



fc2 

2 JO 



d9 0-\{9) + 



= 1 - 0.4409g2 + 0.1316g* 



for g < 1. 



(6.31a) 
(6.31b) 



We see that a'(fc) departs quadratically fro m its a bsolute maximum of unity at the origin. For /c 3> 1 one can replace 
n{9) by its large-0 asymptotic form in Eq. ( |6.30b| ) to obtain 



;(fc) 



1/2 
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40\/26fc2 



for A; > 1. 



(6.31c) 



the 



We see that a;(fc) decays exponentially to zero for large k. For the sake of comparison with the numerical results 
small- and large-A; asymptotic forms for ui{k) are also shown in Fig. ^. 

To summarise, as shown in Sec. VI B the liquid state of a system of randomly crosslinked macromolecules becomes 
unstable when the mean number of crosslinks per macromolecule [M]/N is increased beyond a certain critical value 
Mc/N, corresponding to fi^ = 1. At this critical point the system ex hibits a continuous phase transition from the 
liquid state to the amorphous solid state. As shown in Sees. V^I C| and |VI D , this solid state is characterised by a gel 
fraction q, which grows from a value of zero at the critical point with the classical exponent (3=1 (see Ref. |68| ): 
g~e^/i^ — 1^ {[M] — Mc)/N. The amorphous solid state is further characterised by the statistical distribution of 
localisation lengths 2^~'^p(^~^). In the vicinity of the transition the dependence of this distribution on the control 
parameter e and the (inverse square) localisation length r is determined by a universal scaling function (of a single 
variable) 7t{9), i.e., p{t) = (2/e)7r(2T/e). This universality guarantees that tt{9) need only be computed once for all 
near-critical crosslink densities. As already mentioned in the present subsection, 7t{6) has a single maximum, away 
from which it decays rapidly. Hence, the fraction of localised monomers that are localised on length scales much 
larger than e~^/^ is exceedingly small. Our result for p{t) also predicts that the fraction of localised monomers with 
localisation lengths much smaller than €~^^'^ is also exceedingly small. This provides an a posteriori confirmation of 
the internal consistency of the perturbation expansion in powers of upon which our results rely. However, the 
detailed form of the distribution for localisation lengths much smaller than e~i/2 {e.g., for localisation lengths of the 
order of the radius of gyration of a free macromolecule) is unreliable because such localisation lengths are not within 
the range of validity of the perturbation expansion. The rapid decay of pir) away from its maximum guarantees that 
its moments are finite. This character, together with the scaling form of p(r) ensures that the moments scale in the 
following manner: 

[^-2m] _ {([M]-Mc)/N) . Furthermore, as the distribution has single maximum it is sensible 
to define a typical localisation length ^typ associated with the most probable localisation length. This length ^typ 
obeys the scaling relation ^typ ~ (m^ ~ 1)^^/^. Thus, a simple, reasonable approximation to the true distribution p{t) 
would be a sharp distribution, e.g., S(t — e/2) 



We have seen in Sec. VI D that the order parameter hypothesised in Sec. IIIC and determined in Sees. VIC and 



VI D is a solution of the stationarity condition for the free energy Eq. 
analysed in Ref. [Q, in which it was assumed that q = 1 and p{£,~'^) 



(6.3b) 



This is in contrast with the hypothesis 
2 - ^-^) {i.e., that aU monomers share 



a common localisation length). The hypothesis of Ref. p2| does not satisfy the stationarity condition, and therefore 
only provides a variational bound on the free energy. That our result for the order parameter is a saddle point of 
the free energy, rather than merely a variational bound, is a feature of considera ble sig nificance. The consequent 
exact vanishing of the linear term in the expansion of the effective hamiltonian Eq. ( 5.12bD in powers of the departure 
from the known stationary value streamlines further analysis of, e.g., linear stability, fluctuations, correlations, and 
response to perturbations. 



F. Comparison with numerical simulations 



Extensive numerical simulations of both macromolecular melts and well-crosslinked macromolecular networks have 
been performed by Grest and Kremer, and others; for a review see Ref. On the other hand, until the recent 

work of Barsky and Plischke relatively little attention had been devoted to simulations of the regime in which 
the number of crosslinks is comparable to the number of macromolecules {i.e., the vicinity of the liquid-to-amorphous 
solid transition). It should be remarked that from the computational point of view this is a daunting regime. In their 
simulations, Barsky and Plischke observe a continuous transition from the liquid to the amorphous solid state, and 
extract a universal scaling function describing the distribution of localisation lengths. Whilst their numerical results are 
in strong qualitative agreement with the analytical predictions described in the present article, thus providing support 
for the theoretical picture of amorphous solidification discussed here, there appear to be quantitative differences. At 
present, the precise origin of these differences is unclear: a possible source is the relatively short length of the 
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macromolecules used in the simulations, for which mean-field theory is expected to require substantial fluctuation- 
corrections over a moderately wide range of near-critical crosslink densities. 



VII. INCORPORATION OF DENSITY-SECTOR CORRELATIONS 

Until now we have entirely neglected fluctuations in the fields fl", Qt. and Wk in the computation of the quotient 



of functional integrals in the numerator and denominator of Eq. ( 5.11 ). Instead we have approximated the functional 
integrals by the stationary values of their integrands. This amounts to making a mean-field approximation, in which 
correlations between fluctuations of the fields are neglected. 

As a consequence of this strategy, the excluded-volume interaction has played a subsidiary role: its presence has 
been required, in order to maintain the stability of the physical system with respect to the formation of macroscopically 
inhomogeneous states, e.g., via crystallisation or collapse (which would be detected by their nonzero value of the order 
parameter in the 1-replica sector) , even at crosslink densities large enough to destabilise the liquid state with respect 
to macroscopically translationally invariant states {i.e., equilibrium amorphous solids). However, having accomplished 
this by guaranteeing the stability of the 1-replica sector, the excluded-volume interaction has played no further role. 
Indeed, the precise value of the excluded-volume parameter A does not even feature in the mean-field free energy of 
either the liquid state or the equilibrium amorphous solid state (in the same way that the exchange coupling constant 
of a magnetic system does not feature in the mean-field free energy of the paramagnetic state). Thus, amongst 
the correlations that we have entirely neglected are those between fluctuations of the density of the system, i.e., 
correlations between the 1-replica sector fields fijj . 

There are two reasons why we should seek to improve the theory by incorporating density-density correlations, 
at least at some level. (It should be pointed out that Ball and Edwards have undertaken the task of incorporating 
correlations in their approach to crosslinked macromolecular systems ^,0.) First, at least at high densities, the 
resulting screening of the excluded-volume interaction gives an accurate treatment of interaction effects p^ . Second, 
such correlations would reflect the statistical tendency for the macromolecules to avoid one another. This would have 
the effect of introducing a statistical preference for configurations in which macromolecules stay apart, and a barrier 
between topologically distinct configurations. Thus, one might imagine the following scenario. In a given realisation 
of the supercritically crosslinked system without the correlations associated with self-avoidance there is no statistical 
bias against macromolecules passing through one another. Thus there is no apparent mechanism for creating distinct 
ways for the macromolecules to localise themselves, and the system would exhibit a single family of (translational- and 
rotational-symmetry-related) solid equilibrium states. If, on the other hand, correlations are incorporated, then there 
can be at least statistical barriers between symmetry-unrelated ways for the macromolecules to become localised, 
e.g., topologically distinct interweavings of the macromolecules. Thus, one might anticipate a situation in which 
the incorporation of correlations allows the system to exhibit symmetry-unrelated equilibrium states, a scenario that 
could be revealed through the mechanism of the spontaneous breakdown of the permutation symmetry amongst the 
replicas. Certain ideas in this direction have been explored in Ref. | |33[ |. For a discussion of related issues in the 
context of spin glasses, see Ref. 

A. Free energy 



We now set about the task of computing the quotient of functional integrals, Eq. ( ^.12a ), with improved accuracy 



in order to incorporate the effect of the excluded-volume interaction. With regard to the functional integral in the 
numerator our strategy will be to treat the 1-replica sector in the gaussian approximation, resulting in an improved 
effective hamiltonian for the higher-replica sector field f2j, , which we then treat at the mean-field level. The functional 
integral in the denominator will also be treated in the gaussian approximation. We assume that at the relevant saddle 
points of the impro ved eff ective hamiltonian the order parameter will continue to have the property of being MTI, 



as described in Sec. Ill A 



To implement this strategy, we expand JF„, Eq. ( 5.12b ), to quadratic order in the 1-replica sector field fl^ for an 



arbitrary value of the higher-replica sector field fl^, except for the restriction that ilj, be MTI. Thus we obtain 

ri+l 



a.a'—O k 
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(7.1b) 



in which ly = N/V denotes the number of macromolecules per unit volume. Note that all terms linear in 17^, as well 
as most terms that are quadratic in fig vanish by virtue of the MTI property of fi^. The $7 ^.-dependent correlator 

is discussed further, below, and is computed for a specific form of fi^ in App. Similarly, we expand 
Eq. (5.12c), to quadratic order in the field ojk, obtaining 



5o(|qr 



^^/g-^q.c(^')\)^ 



w 



(7.2a) 
(7.2b) 



The basic Debye correlator 5o(|qP) is compute d in Ap p. ^. Next, we substitute the approximations to Tn and JF'^"", 
quadratic in and Wk respectively, into Eq. ( 5.12a| ), and perform the resulting gaussian integrations o ver and 
ojk- These we do for each value of k in the numerator and denominator, by noting the normalisations Eq. (5.10a) and 
Eq. ( 5.10c ) and applying the well-known result. 



/n n 
J]^ dRe dim exp ( - ^ z"*y^?" 



Q a'— 



/n 71 
J|dRez" dlm2"exp(- ^ z"*y^^"'; 
n „ „ / n 
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Detn A2 
Det„ y^i 



(7.3) 



in which Ai and are arbitrary hermitean complex- valued (n + l) x (n + l) positive-definite matrices and Det^ 
denotes a (replica-space) determinant of an (n + l) x (n+l) matrix. Hence we obtain 



exp ( - ndN^r^i^^^, {0})) w B„ / V^'flcxp 

,2 
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(7.4) 



where /„ is the (n+l) x (n+l) identity on the replica indices. 

Now, our aim is to compute the free energy /, Eq. ( ^.13 ), and to do this we shall need the difference 
$„(/i^, {0}) — $„(0, {0}), as we see from Eq. ( 4.15 ). Thus, in addition to <i>„(/x^, {0}), which we have just obtained at 
the gaussian level of approximation, we shall also need <I>„(0, {0}) at the gaussian level of approximation. This can 

readily be obtained from Eq. (7^) by setting fi^ to zero throughout, except in the measure V £7, Eq. ( 5.9b ), and in the 
coefficient of the quadratic term in the exponent: — Then the logarithmic term vanishes, 0^°" 

becomes (5"" godkp), and A^ becomes A^, and the functional integration can be performed by using the normalisation 
Eq. (5.9b), all /Lt^-dependence cancelling from the result. Thus, we find 



exp I 
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(7.5) 



By forming the difference <i>„(/i^, {0}) — $„(0,{0}) we thus arrive at the following expression for 
[Z({ze, Se] i'e, s'eleli)"] : i'^ which l-rcplica sector density fluctuations are incorporated at the gaussian level, as indi- 
cated by the superscript gdf on T: 
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(7.6b) 
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Following the same line of reasoning made in Sec. VIC, we shall now make a variationa l mea n field approximation 
for the remaining fu nctio nal integral. To do this, we eva luate the effective hamiltonian ( 7.6b ) using the hypothesis 
for Oj. given in Eq. (6.8b). As was the case in Sec. |VI C , this leads to an expression for the effective hamiltonian in 
terms of q and p{t), which we subsequently make stationary with respect to q and p(t). 

The next step in co mput ing [Z (^{ie, Se'yi'e, s'^}eLi)^] concerns the evaluation of the determinant on the replica 
indices, Det„ in Eq. ( 7.6b ). This can readily be accomplished provided we make the assumption that we only 
consider values of that are invariant under the permutation of replica indices (i.e., are replica-symmetric). This is 
indeed the case for the form hypothesised in Sec. Ill C| . With this restriction, takes the form 



which, for a given value of k, has eigenvalues 



1 - (5"" h 



(7.7) 



with degeneracy 1; 
with degeneracy n, 



(7.8a) 
(7.8b) 



so that /„ + X^i^Gk has eigenvalues 

1 + A^ji/(i/k + "^k) with degeneracy 1; 
1 + A^i/(i/k — ^k) with degeneracy n. 

Hence, we find that the determinant in Eq. ( |7.6b ) is given by 



(7.9a) 
(7.9b) 
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(7.10) 



It should be noted that this expression has implicit dependence on n through _ffk, /^k and A^ but not through (7o(|kp). 

Our next step involves making use of the fact that we shall ultimately be taking the replica limit, n — *■ 0. With this 
in mind, we make the following expansions, valid for small n: 



(0) 



(1) 



/^k = /^L"^+ri/i^''+0(n2), 



(7.11a) 
(7.11b) 
(7.11c) 



and note that H^' — go{\k\'^) (as shown in App. By inserting these expansions into Eq. (7.10), taking the 
logarithm, and expanding for small n we obtain 
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Next, we insert Eq. (7.12) into Eq. (7.6b), thus obtaining 
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which is to be made stationary with respect to ftj, of the form ( |6.8b ) on which T^'^^ depends both exphcitly, and also 
impUcitly, through and It should be noted that in obtaining this expression for ^n'^^ {{^j^}) we have made 

use of the result, established in App. that H^^ — go(|kp), without which the limit n — > would not exist. 

We now evaluate for fl^ given by Eq. ( 6.8bD . The second and third contributions to the right side of 

Eq. (k^m) ^^'^ mean-field contributions, which have already been discussed in Sec. VI, and which are computed in 
App. O. The fourth and fifth contributions are fluctuation corrections, resulting fr om th e incorporation of 1-replica 
sector density fluctuations: it is upon these that we now focus. By reorganising Eq. ( 7.13 ) and omitting contributions 
that are independent of both the gel fraction q and the distribution p(t) we obtain 
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where uo 
transition. 



E AgiV/y = y^v. Next, we recall that from Eq. ( |6.11 ) we know that q is small in the vicinity of the 
This, together with the observation that terms for which k = are to be omitted from the fluctuation 
contributions, and that li'^ ^ q (from Eq. (K6) of App. shows that we may expand the last term in Eq. ( 7.14b| ) 
as a power series in liF^ , which yields 



/s'if{g,p} = r"{g,ri 



dvV 



(1) 



oo.go(|k|' 



2di^V 



1 + cro5o(|k|' 



CTo4"^ 



3dvV^ 

k 



l + (7o5o(|k|2) 



Oiq' 



(7.15) 



We now examine the fluctuation corrections in Eq. (7.15) in detail. Equation (K6) shows that h^"^ contains factors 
of the form exp (— fc^/r). For localisation length distributions p{t) that give appreciable weight only to r e, we can 
assert that, where present, such factors effectively eliminate contributions from wave vectors larger in magnitude than 
^/e. Thus, the two terms in Eq. (7.15) that involve h'j^^ respectively scale with e in d dimensions as e^+'^Z^ and g3+'^/2^ 
and thus are negligible for d > 2 and d > 0, respectively. On the other hand, H^'^ contains no such factors. In fact, 
the factorisation of H^^ into t- and k-dependent pieces exhibited in Eq. ( pC5| ) enables us to identify the contribution 
to Eq. (7.15) associated with i?^^' as a renormalisation of the coefficient of one particular term in 
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where the renormalisation parameter £ is given by 
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Now, in App. ^, the functions godkp), 5i(|k|2) ^^^^ g2{\k\'^) are computed. From the asymptotic properties of these 
functions we find that the renormalisation of the coefficient £ is finite for —2 < d < 4, it being infra-red divergent 
for d < —2 and ultra-violet divergent for d > 4. Thus, we find that in the physically relevant three-dimensional case 
the 1-replica-sector gaussian density fiuctuations simply give a finite renormalisation of a coefficient in the variational 
free energy. 

This renormalisation has no effect on the behaviour of g as a function of fM^, which undergoes precisely the transition 
at /z^ = 1 that it undergoes at the mean-field level. Moreover, the only effect of this renormalisation on p(r) as a 
function of fi^, in the vicinity of the transition at = 1, is the finite rescaling: p(r) = (2/e)7r(2T/e) (2/e')7r(2T/e'), 
where e' = e/(l -I- £). However, it must also be recognised that the relationship between the control parameter fj,^ 
and the mean number of crosslinks per macromolecule [M] is modified by the fiuctuation corrections. Thus, once the 
typical localisation length is rescaled, the properties of the system are as found in Sec. VIE. Note, however, that 
although the location of the transition as measured by /z^ is unaltered, the relationship between fi^ and [M]/N is 
changed by the incorporation of fluctuations; therefore, the location of the transition as measured by [M] is altered. 
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VIII. CONCLUDING REMARKS 



As Charles Goodyear discovered in 1839 ||7l|, the introduction of a sufficient density of permanent, random crosshnks 
causes a macromolecular hquid to be transformed into an amorphous sohd, in which a nonzero fraction of macro- 
molecules have spontaneously become localised. This system is a disordered system, in the sense that it comprises 
both thermally equilibrating as well as quenched random variables (i.e., the crosslink locations). In the present article 
we have presented a theoretical description of the physical properties of systems of macromolecules that have been 
randomly crosslinked. Our focus has been on the equilibrium properties of such systems, especially in the regime of 
the vulcanisation transition. 

The qualitative picture of randomly crosslinked macromolecular systems that emerges has the following primary 
features. For sufficiently few crosslinks, the equilibrium state of the system is liquid, thermal fluctuations causing 
all macromolecules to wander throughout the entire container given sufficient time (i.e., all macromolecules are de- 
localised). For sufficiently many crosslinks, the equilibrium state of the system is an amorphous solid state. In this 
state a nonzero fraction of macromolecules are self-consistently localised, exhibiting thermal fluctuations in location 
only over a certain distance scale, which we refer to as a localisation length and which varies randomly from monomer 
to monomer. At a critical density of crosslinks, of order one per macromolecule, there is a continuous thermodynamic 
phase transition from the liquid state to the equilibrium amorphous solid state. The amorphous solid state is most un- 
usual: the mean positions of the localised monomers are homogeneously random, exhibiting no periodicity whatsoever. 
The state is characterised by the fraction of monomers that are localised, together with the statistical distribution of 
their localisation lengths. The symmetry properties of the amorphous solid state are striking. Microscopically, the 
amorphous solid state is one in which translational symmetry is spontaneously broken, certain monomers becoming 
localised about fixed mean positions, in contrast with the liquid state. Macroscopically, however, the amorphous 
solid state retains the translational symmetry of the liquid state, owing to the homogeneous randomness of the mean 
positions of the localised monomers. (Technically, this shows up as the macroscopic translational invariance of the 
value of the order parameter in the amorphous solid state.) This state bears the same relationship to the liquid and 
crystalline states as the spin glass state of certain magnetic systems bears to the paramagnetic and ferromagnetic 
states, in the following sense: the local static density fluctuations in the amorphous solid state correspond to local 
static magnetisation fluctuations in the spin glass state. In both cases these fluctuations vanish, if averaged over 
the entire sample, and hence cannot serve as a global order parameter. Thus, moments of the static fluctuations 
higher than the flrst must be considered. In fully connected long-range spin glasses the distribution of these static 
fluctuations is gaussian (in the replica symmetric theory), so that the second moment — the Edwards- Anderson order 
parameter — characterises the state of the system completely. In contrast, in the amorphous solid state under consid- 
eration here, we find that all moments are equally important. Such non-gaussian statistics are also encou ntered in 
strongly diluted long-range spin glasses |Q , which show a percolation transition as described by Eq. ( |6.11 ) . 

To construct our picture of their physical properties, we have developed a field-theoretic representation of the sta- 
tistical mechanics of randomly crosslinked macromolecular systems. The order parameter capable of distinguishing 
between the various candidate states (liquid, amorphous, crystalline solid and globular) naturally emerges from this 
representation. The presence of quenched as well as annealed variables has been addressed by invoking the replica 
technique. We have derived the saddle-point equation from the effective hamiltonian of this field-theoretic repre- 
sentation, this equation being equivalent to the self-consistent mean-field equation satisfied by the order parameter. 
Whilst it is not apparent how one might obtain the most general solution for the order parameter, we have proposed a 
physically-motivated form for it, which allows for the possibilities of a liquid state and an amorphous solid state. This 
form is parametrised by the fraction of localised monomers, together with the statistical distribution of localisation 
lengths. In fact, this form turns out to yield an exact solution of the saddle-point equation. It should be noted 
that we are only able to proceed with the calculation in the vicinity of the amorphous solidification transition, where 
the typical localisation length is substantially larger than the radius of gyration of an isolated non-self-interacting 
macromolecule. In a refinement of this approach, we have incorporated a certain physically relevant class of corre- 
lations, associated with macromolecular repulsion, at the gaussian level of approximation. No qualitative changes in 
our results stem from this refinement. 

The quantitative picture of randomly crosslinked macromolecular systems that emerges from our field-theoretic 
representation has the following primary elements. At the mean-field level of approximation there is, for any crosslink 
density, a saddle point of the effective hamiltonian that corresponds to the liquid state. However, for crosslink densities 
greater than a certain critical value, this liquid state is unstable. At this critical crosslink density a new saddle point 
of the effective hamiltonian bifurcates continuously from the liquid-state saddle point. This new saddle point, which is 
characterised by a nonzero gel-fraction and a specific distribution of localisation lengths, corresponds to the amorphous 
solid state. The transition between the liquid and amorphous solid states is therefore continuous: in particular, the 
gel fraction and the inverse-square of the typical localisation length both increase from zero linearly with the excess of 
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the crosslink density from its critical value. Moreover, the entire distribution of (inverse-square) localisation lengths 
has a scaling form determined by a universal function of a single variable (which only need be computed once for all 
near-critical crosslink densities). Detailed results for the gel fraction and the distribution of localisation lengths have 
been given. Thus, we see that what emerges is a quantitative picture of the amorphous solid state that confirms the 
qualitative picture proposed roughly a decade ago. 

There are several other contexts in which one can make use of the circle of ideas that we have been using to explore 
the physical properties of randomly crosslinked macromolecular systems. First, one can apply them to a wide class 
of random-network-forming systems. Indeed, a straightforward extension of the present work yields a theory of 
randomly crosslinked manifolds (i.e., higher-dimensional analogues of linear macromolecules |^,Q). Similarly, one 
can address macromolecular networks formed via a random endlinking (rather than crosslinking) process, in which 
one end from each of several randomly selected macromolecules are linked to one other One can also consider 
networks formed via the (freely-jointed) endlinking of rigid or semi- flexible rods jist , subjects that are of particular 
relevance to certain biological structures. 

This circle of ideas has also been used to develop a statistical-mechanical theory of continuous random (atomic or 
molecular) networks, and thus to develop a view of the structural glass transition fl^ . In this case, what emerges 
is a picture of glass-formation in which atomic or molecular units are (permanently chemically) bonded together at 
random, so as to develop an infinite network. Not only do the translational freedoms of the units become localised 
but also do the orientational freedoms. 

A particularly interesting problem arises in the context of NMR data on protein-folding. The basic question is this: 
How much phase space is reduced due to a given number of constraints on the relative positions of certain monomer- 
pairs |7^? We also mention the formal analogy between such questions and Gardner- type calculations in the field of 
neural networks | ]75[ |. Here, again, the basic issue is the reduction of phase space by quenched random constraints, 
even the distribution of randomness being generated in close analogy with the Deam-Edwards distribution, which 
leads to the same formal construction of (n -I- l)-fold replicated systems. 

Finally, we raise the fascinating and difficult issue of the dynamical properties of randomly crosslinked macro- 
molecular systems. In particular, it would be interesting to develop a semi-microscopic theory of the viscosity of the 
liquid state, and its divergence at the onset of solidification, as well as of the dynamics of the solid state itself. Our 
semi-microscopic approach is well-suited to the application of the dynamical lagrangian methods |76| that have been 
developed to address the dynamics of spin glasses p7| . 
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APPENDIX A: REPLICA REPRESENTATION OF DISORDER-AVERAGED OBSERVABLES 



In this appendix we demonstrate how to evaluate the disorder average of products of equ ilibrium expectation values 
by using the replica technique. This complements the demonstration, given in Sec. IV E| , of the connection between 
the disorder-averaged free energy / and the rcplica-Hclmholtz free energy <!>„. 



Consider the disorder average of the product of g equilibrium expectation v alues of physical observables 



[{Oi)^ . . . (Og)x]- (The meaning of the expectation value (• • is discussed in Sec. [HI A| , immediately after Eq. (|3.l|).) 
In terms of sums over disorder realisations and microscopic configurations this quantity is given by 



[{O 



a/xi 



a(x)"Hrie-g^Ai(x)0i---a(x)'Hrge-g.^Ag(x)0g 
a(x)"'tri e-^i Ai(x) • • • a(x)"'trg e-^« A<,(x) 



(Al) 
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Here, tri denotes the trace (i.e., sum) over the configurations one copy (or replica) of the system, without regard for 
distinguishabiUty (i.e., without dividing by any symmetry factor), and tr2 denotes a similar trace over another copy, 
etc. The factor cr(x)(= Ha denotes the exact indistinguishability factor appropriate for the specific realisation of 
the disorder and, accordingly, takes on values between A^! and unity. Ai(x) implements the constraints on the first 
copy of the system, being a factor that is unity for configurations of the first copy that satisfy the constraints x, and 
zero otherwise. A2(x) indicates a similar factor for the second copy, etc. The factor exp ( — Hij indicates the weight 
for configurations of the first copy, and similarly the factor exp ( — H2) indicates the weight for the second copy, etc. 
We now introduce a further (n — g) factors cr(x) '^tr e~^A(x) in both the numerator and denominator, which gives 



[{O 



aixrhri c-^^ Ai(x)Oi • • • ^(x)"'trg e"^' Ag(x)0, 

xa(x)'Hrg+ie-^«+^Ag+i(x)---a(x)~'tr«e-^"A„(x) 

aixyhn e-^^ Ai(x) • • • a(x)~'tr„ e"^" A„(x) 



(A2) 



In principle, this quantity is defined only for integral n > g, for which values of n it is constant. However, in the 
replica approach we regard n as a continuous variable, and formally determine this constant value by computing its 
limit as n — > 0. Thus, we obtain 



[{O 



(a 



9/XJ 



(A3) 



lim [tri e-"'Aiix)Oi • • - trg e"^" Ag(x)Ogtrg+i 0-^'+^ Ag+i(x) • • •tr„e 



A„(x)] 



where we have used the fact that lim„^o(tr e~^A(x))" = 1. As we have seen in Sec. HE, in the context of macro- 
molecular networks the disorder distribution is generated by the partition function itself, which is represented via 
the inclusion of an additional replica. In fact, it turns out that the naturally emerging order parameter involves the 
expectation value of variables associated with this additional replica, as we see in Sec. [II. Such expectation values 
may be regarded as introducing additional x-dependent factors into our disorder averages, e.g., 



[e(x)(Oi>x...(Og>x] = 

lim [e(x)trie-^iAi(x)Oi---tr<,e-^'Ag(x)0<,tr<,+ie-^'+iA<,+i(x)---tr„e-^"A„(x)] 



(A4) 



By using Eq. ( |A4|) and th e Dea m-Edwards distribution, Eq. ( 2.12 ), in a manner analogous to that used in the context 
of the free energy in Sec. [VB, we arrive at the result 



[e(x)(o 



lim ■ 



Oo{cUs)}0,{cl{s)}---0,{c!{s)} 



Q = 



E 

n+1 

E 

n+1 



(A5) 



where (Oo)x = ©(x)- 

We now focus on a particularly important example of a disorder-averaged observable, viz, the amorphous solid 
order parameter Eq. (^^). This quantity is obtained as the special case of Eq. ( A5) for the choices: 



e(x) = (exp(ik0.c,(s))), 
Oi = exp (ik^ • Ci(s)), 

Og = exp(ik9 •c,(s)). 



(A6) 



By using Eq. ( |A5| ) for this case we see that 
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1 f 

Y X! / ^^^^ ■ ^^(*)))x(exp (ik^ • c,(s)))^ • • • (exp {ik^ ■ c,{s))) 
i=i •'^ 

lim ^-^^ — ^ , (A7a) 

-.'(^ E n 



where Q^. is defined in Eq. (4^), and I — {k", k^, . . . , k^, 0, . . . , 0}. As was anticipated in Sec. [II C| , the order 



parameter probes not only the g physical copies of the system but also the additional copy used to generate the 
crosslink distribution. 

APPENDIX B: WIENER CORRELATOR 

In this appendix we derive the basic correlator associated with the Wiener measure: 

r ^ 1 ^ 

(|exp(-i^kp-c(sp))^^ = (jj^-^,^ exp ( - - ^ Spp-kp-kp,y (Bl) 

where Sppi is a function of the pair of arclength coordinates Sp and Spi defined via 

Spp> = min (sp,Sp/) . (B2) 
In terms of the Wiener measure, the correlator is given, up to normalisation, by 

exp ( - i^kp • c(sp))^ oc /x>c exp (^-i / ds|c(s)|^) cxp (^-i ^ k^ • c(sp)) , (B3) 
p=l ^ •' -^0 p=l 

where the overdot denotes a derivative with respect to s. We express the configuration of the macromolecule in terms 
of the position of the end at s = 0, together with the tangent vector c(s) via: c(s) = c(0) + ds' c{s')9{s — s'), where 
0{s) is the Heaviside f?- function. Then the measure "Dc is given by l?cdc(0), and the correlator becomes proportional 
to 

1 1 
/ dc{0) exp ^— zc(so) '^kp^ / Vc exp ^— - / ds\c{s)\^^ exp 

P=l P=^ 

By performing the integral over c(0) we obtain the Kronecker (5-function factor: ■ By performing the 

integrals over the tangent vectors c(s) we obtain the gaussian factor: 

(B5) 



exp ~ 2 $Z ' ^P' / '^^ ds'9{sp — s)0{spi — s')S{s — s')^ . 
p,p'=i ° 



By performing the arclength integrals, and by setting to zero the wave vec tors {kp}p^]^ in order to establish the 



correct normalisation factor, we obtain the Wiener measure correlator Eq. ( |B1[) . It should be noted that because 
min (sp, Sp') = ^{sp + Spi) — ^\sp — Sp'\, and because of the Kronecker ^-function factor, the exponent of the Wiener 
measure correlator can also be expressed as p,^-^ kp ■ kp>\sp — Spi\/4. 

APPENDIX C: DEBYE FUNCTION AND RELATED FUNCTIONS 

In this appendix we give the basic Debye function godkp), as well as the two related functions gi(|kp) and .g2(|kp). 
The Debye function go(|kp) is defined as the integral over arclength variables of the Wiener correlator Eq. ( |Bl[ ) for 
the case r = 2 and — ki = k2 = k: 
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9oi\ 



kf ) = dsids2 (e^*" i^'^^-'^i-^)) = dsids2 e-''' 



-fc2|si-S2l/2 



r2/45, ifP^l; 

\32/3fc6 ifP^l. 



(CI) 



The function f;i(|kp) is defined as 

gidkp) ^ £ d.id.2d.3 e-'=^l---l/^( - 5i,3 + 52^3)' - { ;[ J I; (C2) 

The function (72(|kp) is defined as 

JQ 



(C3) 



fn d = 3 the relevant correction due to gaussian density fluctuations involves the quantity gi(|kp) — (72(|kp), the 
asymptotic behaviour of which is given by 

»,(wv».(ikP)~{;/;§i, '[^^j;: (C4) 



APPENDIX D: EFFECTIVE HAMILTONIAN EVALUATED FOR THE ORDER PARAMETER 

HYPOTHESIS 

1. Quadratic contribution 



We now compute the contribution in the n — » l imit to the effective hamiltonian, Eq. (5.f2b), that is quadratic 
in ri^, for the specific form of flj^ given in Eq. (3.14). No approximations will be made. To this end we focus on the 
quantity 



f J ;F^t 



hm -^Y. I^l'- 



(Dl) 



Inserting flj^ from Eq. (3.14), we obtain 



yn 



= q' 



oo 

J J dTip{Ti) j dT2p{T2)5{T - {t^^ ^ T2^) 



(D2a) 



We now add and subtract the terms in the 0- and I-replica sectors in order to relax the constraint on the summation 
over k indicated by the overbar on the summation symbol. In fact, owing to the factor of ^''p, the summand vanishes 

for k in the I-replica sector, so that we obtain 
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V 



-^ + q^ J drj dnpin) J dT2p{T2) s(t - (rr^ + r^') ') 

x^E4>-p(-^V2r). 



(D2b) 



We now focus on the remaining, unconstrained summation on the right hand side of Eq. ( D2b ) , which we compute 



via the following steps. First, we introduce an integral representation of viz. 



(D3) 



where the integral is taken over the volume V. Then we convert the summation over k into an integral by using 



yn+i J 



dk - • ■ . 



(D4) 



where = 00=0'^'^" andcfk = (27r) "^dk, and the integral is taken over the entire k space. Thus we obtain 
V^. E 4'i ( - ^V2r) ^Jdke^p[~ P/2r) dc exp (^c • J^^^ k^J 



/^.c{/.?kexp(-k.k/2r)exp(.k.c) 
J dcexp ^ — (n + l)rc • c/2 



(n+l)d/2 



(r/27^)"''/'(l + n)-'^/^ 



(D5) 



where we have used the gaussian integral and, in the last step, have assumed tha t r <C F^/'^ for inverse square 
localisation lengths that are given significant weight by the distribution p{t), Eq. (3.14). Thus we find 



<?■ , 2 



dr dTlp{Tl) (iT2p(T2) (5(t - (T]"^ + T2"^'' 



x(T/27r)"'/'(l + n)-'^/2. 



(D6) 



As we shall need this quadratic term only in the vicinity of n = 0, we expand for small n using Eq. (4.14a), and by 
taking the n — > limit we obtain 



1 9 1 

lim 



pea POO fy'^/d _i 

' dn p{Ti) dT2 p(t2) In (rf 1 + T^^) 



(D7) 



2. Logarithmic contribution 



We now compute the contribution in the n ^ limit to the effective hamiltonian, Eq. ( ^.12b| ), that can be ide ntifie d 



and the partition function of a single macromolecule coupled to fi^ for the specific form of fi^ given in Eq. ( 3.14 ) 
The calculation will be undertaken as a perturbation expansion in the typical inverse square localisation length, to 
first order in this quantity. Thus, we focus on the quantity 

exp(^2^V-"^.Ref]|y dse'''-^^'^jj^^. (D8) 



First, we observe that Qj^ has only been introduced for k ■ h > 0, as follows from the discussion after Eq. (5J). We 
are therefore free to introduce Qj^ for fc • n < at our convenience, and we do so via the definition 
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=n* r for fc • n < 0. 



(D9) 



This allows us to eliminate the Re operation and to extend the range of the summation in Eq. ( [D§| ), which becomes 



exp / ^'^ 



W 



n+1 



Next, we insert VLf_ from Eq. (3.14), which gives 



-ifc-c(s) 



W 



n+1 



(DIO) 



(Dll) 



We now add and subtract the terms in the 0- and 1-replica sectors to the summation over k in order to relax the 
constraint indicated by the overbar on the summation symbol. In fact, owing to the factor of '5^''^, the summand 

vanishes for k in the 1-replica sector, so that we obtain 



exp ( - ^i\V'") ( exp (^i\v-^^ €i r '^^^^^^ "^"''^'^ /' ' 

I' "0 



-ik-c{s) 



W 



n+1 



(D12) 



The next step is to make the power series expansion of the exponential in the expectation value, and make r-fold use 
of the integral representation of the Kronecker ^-function, 



dc exp ( ic • k J , 



(D13) 



in which the c-integral is taken over the volume V, to obtain 

t 



exp U^qV^'') ( exp f 2^V""^ Re^l / dse 



y ynr j.i 

''=1 ' ku---,kr- 



w 



n+1 



cxp^i 2_^Cp l^kp 

r -, n 



/ dTip{Ti) ■ ■ ■ dTrp{Tr) CXp ( - - ^ — ^ | k^ | 

«1 r n 

X / dsy ■ ■dsr(exp ( — kg ■ c°( 



W 



n+1 



(D14) 



The remaining expectation value in Eq. (D14), ( exp ( — i X]p=i X]a=o ' '^"('*p)))I^i' factorises on the replica index, 

giving ria^o ( ^^P ( ~ ^J^p^i^p ■ '^('5p)))i^- Each factor in this product is of the form of the expectation value 
computed in App. in terms of which we express the remaining expectation value in Eq. (D14). The result contains 
an [n + l)-fold product of Kronecker (5-functions, which we replace with the integral representation 



TT s^''^ - - 



„ n r 

— dm0...dm"exp(-z^m".^kg 

n 1 



(D15) 



in which each of the (n+1) m-integrals is taken over the volume V. Next, we convert the summations over \ ki , ■ • • , ^r} 
to integrations by using Eq. (PJ) r times, after which the summation on the right hand side of Eq. (D14) becomes 



» n r ^ r n 

/<ffci...<?fc.exp(*^^kg.(cp-m")) cxp(-- ^ 7^2 ^ kg • kg, 

a—Op—1 p,p' — l Q— 



1 f dci dCr 



dxnP ■ ■ ■ dm" 
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J dm J dki ■ ■ -dkr exp (^i • (cp — m) 



X exp 



1 



p=i 



T?^''^ k • k , 



p,p'=i 



n+l 



(D16) 



Here, Ti-^pl is an (r x r)-matrix-valued function of the r arclength coordinates {szyj^^i and the r inverse square 



pp' 



locahsation lengths {ti,}^,^^, defined in Eq. (Ff ). We now focus on the quantity in this expression that is raised to 
the {n + power. 

By performing the gaussian integrations, first over {ki, . . . , k^} and then over m, we obtain 



J dm J dki ■ ■ ■ dkr exp (j. k^ 



[Cp - ml exp 



1 



'^pp' ^P ^P 



(27r)"('^-i''^/2 (w*"^) Det^'') n'^< 



p,p'=i 

-d/2 

exp 



9 pp' 



P,p'=i 



(D17) 



where we have introduced the (r x r)-niatrix-valued function C^"^^ of the r a rclength coordinates {siy}^^^ and the r 
inverse square localisation lengths {r^jj^^i, which is built from TZ^^\ Eq. (Fl), in a manner described in App. The 
gaussian integrals in Eq. ( D17| ) are convergent, o wing to the positive-d efinit eness of the eigenvalue spectrum of T?.'-''^ 
and of yV'('") for finite t. By inserting the result (D17) into expression (D16) we obtain 



"^X!^~7r~_^ dsl---dSr J dTip{Ti)- ■ - dTrPiTr)^ J 



1 f dci dCr 



(2^)-(-l)("+l)rf/2 {y^ir) ^(.)^ 



n + l 



(D18) 



The quantity W^'') is built from Eq. (Fl), in a manner described in App. ^ 

Next we perform the integration over {ci, . . . , c, }. This integration is not quite gaussian, instead being quasi- 
gaussian, owing to the presence of a single zero-mode, the eigenvalue spectrum of C'-'"-' containing a single zero 
eigenvalue, the remaining r — 1 eigenvalues being positive-definite. The presence of this zero-mode can readily be 
ascertained by observing that from the definition of C''^\ Eq. (F2c), we identically have X]p'=i^pp' — 0, i-e., the 
normalised r- vector (1, 1, . . . , 1)/y^ is an eigenvector of C^^J, with zero eigenvalue. The necessary quasi-gaussian 
integral is computed in App. ^. By using it, Eq. (D18) becomes 



OO 9^ r 



-y 

r— 1 



dsr ■ -dsr / dTip{Ti) ■ ■ ■ dT^p(T^)(27r)-('^-i)"''/2 



X (W('-) Det^--) 7^('■)) {n + ly^^-'W^ ^^-ly^ir) j^^^M ^(.) ^^^C^) ^(r)^^ 



-d/2 



(D19) 



~ (r) 

where Det C^^-* indicates the quasi- determinant of C^^\ i.e., the product of all the nonzero eigenvalues of C^''-'; 
see [Q. We now make use of the result, established in App. that the factor W^'") Det^'') 7e('-) Det^''^ C^'') IS 
identically unity. Thus, Eq. (D19) is simplified, becoming 



— — / dsi - ■ -dSr / dTip{Ti) ■ ■ ■ dTrp{Tr) 



(D20) 
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We take this expression and insert it into Eq. (D14) to 

obtain an expression for exp ^2/i^l/~"^|,Re fii ds e**^ '^*^^-'^ ^ that is vahd for n > —1. By expanding this 
result for smaU n, we finaUy obtain the desired quantity: 

•c(s) 



Urn ^ln(exp (2/^2^~ny-^ j^g^* /■ ^^^fk- 

= (c-^'« - (1 - At'g)) hi (^2/d/2^e) - e-'^'« ^ 



r— 1 



X ^ dsr ■ -dsr dTipin) ■ ■ ■ dTrp{Tr) hi (w^'^^ Det^''^ 7^('^)) . 



(D21) 



APPENDIX E: QUASI-GAUSSIAN INTEGRATION 

Consider the foUowing integral, 

" ' ' "cp.Jp), (El) 



^ / dci • • • dc^ exp ( - i ^ C^^! Cp-Cp,^exp(^-iJ2' 



taken over r copies of the volume V , in which C^^-* has as its sole non-positive definite eigenvalue the zero eigenvalue 
corresponding to the r-eigenvector (1,1,..., l)/\/r. We shall need both the general case, in which the sources {ipYp=i 
are arbitrary, and also the special case, in which the sources {Jp}p=i all vanish. By working in a basis in which C^^^ 
is diagonal and assuming that V is sufficiently large (or, equivalently that no nonzero eigenvalue of C^^-* is arbitrarily 
small) one finds that this quasi-gaussian integral is given by 



where C*-''-' is the quasi-inverse of C*-''-', i.e., the eigenvector expansion of the inverse of C^'-* from which the term 



corresponding to the zero eigenvalue has been omitted; see |^^. Similarly, as mentioned in Sec. D2, Det C*^''^ 
indicates the quasi- determinant of C'-'"'; see [Q. The factor of r'^/^ is subtle, but is familiar from the context of 
so-called collective coordinates. It arises from the fact that owing to the presence of the zero eigenvalue there are d 
integration directions for which convergence is not controlled by a gaussian integrand. Not only do the corresponding 
integrations yield a volume factor; they also each yield a factor of r^/^ by virtue of the limits on them determined by 
the form of the corresponding eigenvector. For the special case in which the sources {Jp}p=i all vanish we have: 

1 |^dc,...dc.exp(~i ± C^cp-v) .r''/^(2.)(^-^^^/^(Det^'-^CW)''''. (E3) 
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APPENDIX F: PERTURBATION EXPANSION AT LONG LOCALISATION LENGTHS: FREE ENERGY 

Consider the quantity 7?,'''^ an (r x r)-matrix-valued function of the r arclength coordinates {s,y}J^^i and the r 
inverse square localisation lengths {t^}^=i, given by 

K1-^^p'^pp'+Spp>, (Fl) 



where S is given by Eq. (B2). We have found it useful to construct from TZ'^'^^ several auxiliary quantities: 
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- V 

p'=l 




7 

pp' 




r 

p=i 






Mr) 









(F2a) 

(F2b) 
(F2c) 



We shall need to develop perturbative expansions in {t^Yu=i of ^'^ (yV^^^ Det^'^-' TZ^^^j to linear order and of C*^''^ to 
quadratic order. To this end, we introduce the (r x r) identity matrix X^^^ = Spp' , and make the definition 



7^, 



(r) 



pp' 



'p "pp ' 



so that (7?., 



(r)N-l 



= TpSppi. First, we consider InDet^'"-' T?,^'') 
InDet^'') 7^('^' = InDet^'') (7^^''^ + 5 

= lnDetW7^^'■^ (xW + (7^^'^)) -5 
= InDet^"^) 4'') + InDet(^) (x^^) + (4''^)"' • 5) 
= InDet^"^) 7e[,'') + Tr^'') In (x^'') + (7^^'^)) • S) 

where Tr'-''-* denotes the trace of an rxr matrix. Second, we consider {vS"^^^ ^ . From Eq. (Fl) we have TiP^ 
from which follows the Dyson-type equation 

Iterating the Dyson-type equation once and then truncating gives 



pp 



TpSpp' ~ TpSpp'Tp' -f 0{t^). 



By using this result in Eqs. (F2a) and (F2b) we obtain 

r 

U^p^ ^Tp-Tp^Spp^Tp,, 



p' = l 



p=i p,p'=i 



Third, we use Eqs. ( |F4| ) and ( [F7bD to obtain 

Sp=i "^p 



In 



Det*'') 7^('': 



= In 



n;=i 



p=i ''p 



E '^p^pp 
\p=i 



Sp,p'=i ''p ^pp' ''p' 



Finally, by using Eqs. ( |F7a| ) and ( |F7bD in Eq. ( |F2cD we obtain 

TpTp, 



= ( TpSpp. 



^p ^ pp' '^p' '^P^f^ 



pv I 



+ 0(r3). 



(F3) 



(F4) 



(F5) 
(F6) 

(F7a) 
(F7b) 

(F8) 



(F9) 



48 



In reducing the free energy ( ^.1C| ) to the form (6.14) we have used, inter alia, the perturbatfon expansfon (F8). In 
this way, the free energy reduces to an assembly of terms each being a functional of p{t). Each term has a coefficient 
determined by integration over the arclength variables {si, . . . , Sr} of integrands arising from factors of Sppi , defined 
in Eq. (B2), which depend on the arclength variables. The necessary integrals are readily expressed in terms of the 
following ones: 



/ dsiSii= / dsi min (si, si) = / dsiSi^l/2, 
Jo Jo Jo 

pi pi pi pi pi psi 

/ dsi / ds2Si2 = / dsi / ds2 min(si,S2) = 2 / dsi / ds2 S2 ~ 1/3. 
Jo Jq Jo Jo Jo Jo 



(FlOa) 

(FlOb) 



APPENDIX G: A USEFUL IDENTITY 



We make repeated use of the identity 

W^"") Det^'') 7^('^) Det^''^ C'"") 



(Gl) 



where TZ''^^\ yv('") and C'^''-' are, respectively, defined in Eqs. (|F1|), (F2b) and (F2c) of App. ^. To derive this identity, 
we evaluate the following quantity [which arises, e.g., in Eq. (D16)] in two ways: 



^ dm dki- -dkr exp (j,^kp ■ (cp - m)^ exp - i ^ "^pp'kp • kp- 
1/ p=i p,p'=i 



(G2) 



First, by integrating over {cpjp^j^, then over {kp}p^]^, and then over m we obtain the result: unity. Second, by 
integrating over {kp}^^j^,then over m, and then, by using the quasi-gaussian integral Eq. (E3), over {cpjp^j^, we 
obtain the result: 



V-lwWDetW7^MDet^'■^CW 



~d/2 



(03) 



hence the identity Eq. (p|). 



APPENDIX H: LAPLACE REPRESENTATION OF FREE ENERGY 



In this appendix we describe in detail how to exchange the dependence of the three contributions to /™ , Eq. (3.14), 
from p(r) to its Laplace transform p(f). We begin by noting two integrals: 

'^^ ' -'^) , (Hla) 







T 



1 



dfe-^^. (Hlb) 



The latter integral is elementary; the former is an example of a FruUanian integral; see Ref. ||79|]. We use Eq. ( Hla ) 
to express the first contribution to Z™"^ as 

= y dnpin) J dT2p{T2) j ^ (-e-^ - e-^(^i+"^) + e-^"i + e-^^^) 







dr 

— {-p{ff + 2p{f) - e--) , (H2a) 
r 



where the curly braces {• • -j^ were defined immediately following Eq. (6.14). By following the identical strategy, we 



use Eq. ( Hla ) to express the second contribution to /™'' as 
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{in (^^^^) } = {In in +T2+ T3) - In (n) - In (ra) - In (ra)}, 

— {-p{ff + 3]5(f) - 2e-^) . (H2b) 
To express the third contribution to Z™"^ in terms of p{t) we make use of Eq. ( plbl ). This yields 

\ ^, ^ ^ / C?-rip(Tl) / dT2p{T2) \ ^ 

nOQ nOQ nOO pOQ 

= / dnp{Tl)Tl / dT2p(T2)T2 / dfe-^^^'+^-'^ = / df{dp/dTf. (H2c) 

7o Jo Jo Jo 

We now take the functional derivative with respect to p(f) of these three contributions to y™"^. Being local, the 
first two are straightforward to compute, respectively yielding 



In 



5p(t) \ \ T1T2 



Ti + r2 M 2 



T 



(l-p(f)), (H3a) 



' ^nflL±^H ^^fl-p(ff). (H3b) 



SpiT") I V T1T2T3 

To evaluate the functional derivative with respect to p{t) of the third contribution to /™'' requires an integration by 
parts. The integrated piece vanishes because p{f) is to be varied at neither f = nor f = 00, due to the boundary 
conditions. Thus we obtain for the third contribution to J™"^: 

1 (H3C) 



Sp{f) [ti +T2 J ^ C?f2 



APPENDIX I: ORDER-PARAMETER WEIGHTED AVERAGES 

In this appendix we focus on the computation of the following quantity, defined for arbitrary I and V: 



exp i^^i^qV-'' I] y e-^'=-^(") 6*^1 J drpir) exp ( - P/2t] 



w 

n+l 



(11) 



In addition, we make two applications of it. We proceed by expanding the exponential, which yields 







+1 ^1 V^'^-rl Jo JO 



p=l 



n+l 



(12) 



Next, we observe the factorisation of the Wiener measure correlators on the replica index, use the expl icit res ult f or 
this correlator given in App. |b[ use integral repre sent ations for the Kronecker (5- functions [Eqs. ( D15 ) and ( D13 )], 
and convert summations to integrals by using Eq. (D4), thus obtaining 



^("''t'^) f\t C dt' e-\'-''\^ 
'+'''0 io io 

°° ,,2r„r /•! 



-—— 2_. j — / dsi- ■ ■dSr+2 / dTip{Ti) ■ ■ ■ dTrp{Tr)— / dci . . . dCr / dm" • • • dm" 

"I"- Jo Jo ^ Jv Jv 
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/r n n r r 

dh ■ ■ -dkr exp [^Y.''p■Y.K) (- * E (1" + 1'" + E K)) (-E 

- n r 71 r n r 

^ ^^P (- 2 E E -k^' cxp (- E E K Sr+i,p) exp (- E E 1'" ' 5, 



■r+2,p 



a— p,p' — 1 



X exp — -/ S; 



?r+l,r+l 



Q=0 p=l 

exp ( - -f'^5,-+2,r+2 ) cxp ( - r Z'5,.+ i,r.H 



a=0 p=l 



(13) 



We now recognise that the integrals over {m"}J^^g and {k", . . . ,k"}2^Q factorise on the repHca index, giving 



l+l'fi 



/•I . 1 °° ,,2r„r pi 



"'0 



^ — t J- 

r—1 



dsy ■ ■dSr+2 I dTip{Ti) ■ ■ ■ dTrp{Tr 

Jo 

1 f 

— / dci ■ ■ ■ dCr 



X exp ^ — -P Sr+l,r+l — 2^'^ '5r+2,r+2 — 1 ' 1' Sr+l,r+2 

^ n ( / dm dki - -dkr exp (^i'^Cp- kp^ exp - im • (r + 1'" + ^ kp 
Q=o ^-^^ p=i p=i 

X exp ( - ^ k^/2rp) exp ( - - ^ kp • kp, 5pp' 



p=i 



p,p'=i 



r r 

exp ( - ^ P • kp Sr+i,p) exp ( - ^ 1'" • kp S, 



■r+2,p 



p=l 



p=l 



(14) 



At this stage we focus on the factorised integrals over {k"}^ and m" occurring under the product. As they are 
gaussian integrals they can straightforwardly be performed (we prefer to integrate over {k"}p^]^ first and m" second), 
yielding 



(2.) 



-(r-l)!i/2 



X exp 



(W^*^) det7^('■)) "^^^xpl^- |P + 1'"|V2W^''^ 

~ 9 E (^P + + '^'•+2.P) • + '5.+l,p' + ll'" Sr+2,p') 



P,P' = i 



X exp 



(P + 1'") • ^ (cp + zP 5,+i,p + 5,+2,p)) , 
p=i 



(15) 



where U^''\ W^'') and C^''^ are respectively defined in Eqs. ( |F2a| ), (|F2b| ) and ( |F2q ). Next, we insert the result of 
integrating over {k"}^^^ and m" into Eq. ( ^ and focus on the remaining integrals over {cpYp^i. Just as we 
encountered in App. y, this integration is not quite gaussian, instead being quasi-gaussian, owing to the presence a 
single zero in the eigenvalue spectrum of C^''-'. The necessary quasi-gaussian integral is computed in App. By using 
it, Eq. (U) becomes 

f'dt /'di'e-l*-*'l'" 



^0 



+4+I',0"f7^E^~7r~ / dsvdSr+2 J dTip{Ti) ■ ■ ■ dTrp{Tr) 

X (27r)-^'''-'^"'/'(w('-) DetW 7^W)-"'/'(n + i)-^'-!)'^/^ 
X (r~i Wt"^) Det('-) U^^^ Det^''^ c^r)yd/2 

X exp ( {W'^-y^Y.K''^ ^r+l,p + i'^Sr+2,p + I ■ i' iSr+l,p + 5,+2,p) 
^ p=l 

^^p (2^' E CpJ' ^r+l,p Sr+l,p' + -i'^ ^ CpJ, Sr+2.,p Sr+2,p' + I ' E '"pp' ^r+l,p Sr+2.,p' 



X exp 



p,p'=i 



p,p'=^ 



p,p'=i 
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X exp 



(- (z2 + f'2 + 2f-/')/2W^''^) 

J, 

exp ( - (n + 43(l'5r-+i,p +i''5.+2,p) • (i5,+i,p. +i'5,+2,p')/2 



p,p'=i 



(16) 



We have simplified this result by making use of the identity given in Eq. ( |Gl| ) of App. O. We have further simplified 
it by observing that as C*^''^ is the quasi-inverse ||7^ of C'-''^ (see App. ^ and footnote ||78|| ) , the relevant zero-mode 
being (1, 1, ... , 1)/^^, we have 



pp 



(I7a) 
(I7b) 



In addition, we have used the fact that there is a single zero- mode to ascertain that Det^ ■* ((n -I- 1) C*^''-') = (n -I- 
l)'^-iDet^''^CW. 

As our first application of Eq. (|^), we set I ^ I' = 0, thus obtaining a normalisation factor that we shall use 
subsequently: 



^ + T7^y^. i~ / dsi---dSr+2 dTip{Ti) ■ ■ ■ dTrp{Tr) 

yn ^ ^! 

X (2^)"<''"'^"'/'(w('-^ Det« 7^M)""" 



1) 



-{r-l)d/2 



(18) 



As our second application of Eq. (|^) we compute the right hand side of Eq. ( |6.26| ) by forming the quotient of E q. (|6 ) 
with V — 6 and Eq. (|8|). By making use of the indentity Eq. (Gl) and taking the limit n — > we obtain Eqs. (6.27a) 
and ( l6.27b[ ). 



APPENDIX J: PERTURBATION EXPANSION AT LONG LOCALISATION LENGTHS: ORDER 

PARAMETER 

In th is ap pendix we compute the per turb a tive e x pansi ons of T^"'"-^ and T^^^ needed to compute the right hand side 
of Eq. ( |6.29| ). By using the definitions (^), {^2^, ( p2b| ) and (^2c|), we find 



T(2) =Ti+T2+0(t2). 



(Jla) 
(Jib) 



By inserting these results into Eq. (6.29), and using Eq. (6.28), we obtain 

qp{T) ^ {I- q)n^q j dsids2 dTip{Ti) 5{t ~ ti + t^\si - 82]) 
Jo Jo 

1 f°° 

+ / dsids2ds3 / dTip{Ti) dT2p{T2)S{T - Ti - T2) + 0{e^ 

^ Jo Jo 



(J2) 



Next, we expand the Dirac (5-function, (5(t — ri + t^|si — S2I) « (5(t — ri) + t^|si — S2\S' (t — ri), and perform the t 
and s integrals. (Equivalently, we take the Laplace transform of Eq. (p^, expand the resulting exponential function, 
perform the r and s integrals, and back-transform the resulting nonlinear ordinary differential equation.) Finally, we 
transform to the scaling form 7r{6) defined in Eq. (6.21a), and observe that it satisfies Eq. ( |6.25a ). 



APPENDIX K: CORRELATOR FOR ONE-REPLICA SECTOR FLUCTUATIONS 



In this appendix we compute the quantity Qq defined in Eq. (7.1b) for Qj^ taking the form given in Eq. (B.14). 
Ultimately, we shall be concerned with the behaviour at small n of the diagonal and off-diagonal elements of via 



the quantities H^^ and h'j^^ defined in Eqs. (7.11a) and (7.11b) 
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First, we note that by using the invariance under k — > ~k of the hypothesis for fij,, Eq. (3.14), along with the 

adding and subtracting of terms in the 0-rephca sector, so as to relax the constraints on the summations in the 
exponents in the numerator and denominator of Eq. (7.11), we arrive at the form: 



Jo ^ 

X exp [f^^qV-^'Y^ dse-'''-^'-''^d'^l drpir) exp ( - fcV2r)) ^ 



w 

n+l 



e^pL^qV-^'Y, rfse-*'^(")4''W dTp(T)exp(-P/2T 



w 

n+l 



(Kl) 



Second, we set I = {— q, 0, . . . , 0} and I' = {q, 0, . . . , 0} to obtain the repHca-diagonal elements of the numerator 
correlator in Eq. (Kl): 



5o(|qP) + -^X! dSi---dSr+2j dTip{Ti) ■ ■ ■ dTrp{Tr) 

X (27r)-^'^-'^"'/'(wW DctW 7^M)-"'/'(n + ly^^-'^"'^ 



X exp 



1 + Sr+2,r+2 ^ 2Sr + l,r+2) 



P,P' = 1 



(K2) 



Third, we set I = {— q , 0, 0, . . . , 0} and I' = {0, q, 0, . . . , 0} to obtain the replica-ofF-diagonal elements of the numerator 
correlator in Eq. (Kl): 



Sql + Tl^Y] 1— / dsi-- ■dSr+2 / dTip{Ti) ■ ■ ■ dTrp{Tr) 

X (2^)"('^-')"''/'(wW DetM 7^W)""'/'(r^ + ly^'^-'^"^' 



exp ( - 7;Q^{{Sr+i, 



r + l + Sr+2,r+2) 



p=l 



— ^ Cpp' {Sr+l,pSr+l,p' + Sr+2,pSr+2,p') ^ ^ _^ ^pp''5'-+l,p'5r+2,p 



P,P' = 1 



P.P' = 1 



(K3) 



and off-diagonal elements, H^^^ and /ij^^'', of Q^, Eq. (Kl). At this stage, all dependence on n is explicit. Thus, by 
making series expansions in n about n = of the denominator and respective numerators we are able to identify H^'^ , 
H^^^ and h'^'^ of Eqs. ( |7lla| ) and ( [7llb| ) (for q 7^ 0): 



,(0) 



Next, by using the denominator, Eq. (|8|), and respective numerators, Eqs. (K2) and (K3), we build the diagonal 



<^=5o(|qn, 

H'^^ = e"^ ^ / , exp ^ - q^\sr+l - Sr+2|/2^ / dTip{Ti) ■ ■ ■ dTrp{Tr) 



^(9^/2) ^ C^p' ( — '5r+l,p + 5r+2,p) ( — 5r+l,p' + 5,.+2, 
P:P' = 1 

2r r /•! /"Oo 
/,(0) = e-p^9 ^ / rf^^. . .^^^^2 / . . . dTrp{Tr) 

r=l ^' "'0 Jo 



P'J' 



(K4a) 



(K4b) 



/I 1 1 

exp f - g2 (5^+1^,.+! + 5^+2,r+2) + Y^ - ^p'^ {^r+l,p + S, 

^ P=l 



■r+2,pj 
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- E C^pl'Sr+l,pSr+2,p')], (K4C) 

where we have restricted our attention to q 7^ 0, which is aU that is necessary, and have used the fact that C^^^ 
vanishes identically. 

As a final component of this appendix, we expand iJq^^ and /iq°\ assuming that the gel fraction q and the inverse 
square localisation lengths to which p(t) gives appreciable weight are both asymptotically of order 1), i.e., which 
is appropriate in the vicinity on the transition, as we can verify a 'posteriori. The accuracy t o whi ch ifq^' and 
must be computed can be established by observing the manner in which they appear in Eq. ( 7.15| ), from which we 

see that can be computed to O^e") by computing to O(e^) and /iq^-* to C(e^). 

We first focus on i/q^"*. As we see from Eq. (^), C^^! is of order r, so that we need retain only the r = 2 term in 
Eq. (K4b). Thus, to sufficient accuracy we obtain 

^q'^ = T?'/?' / ds3ds4exp ( - g2|s3 _ S4I/2) 

4 Jo 

X / dsi_dS2 ( " Sz.p + 54,p) ( - 53,p' + S^^pi) \ TpSpp' ~ ^r'^'' — > 

Jo p^p,^^ L Z..=i^aJ 



2mV|;^^ 9^(.9i(|qn-52(|qn), (K5) 



1 4 2/ '^l'^2 



where the curly braces {• • -j^ were defined immediately following Eq. ( 6.14 ), and the functions ffidqp) and 52(|qP) 
are given by elementary integrals defined in App. |^. Note the convenient factorisation into q-dependent and p(t)- 
dependent terms. 

Last, we focus on /iq""*, for which we also need only retain terms with r < 2. However, the r = 1 term is not 
identically zero, and as it carries only a single power of q care must be taken to compute it with appropriate accuracy. 
(In practice, the absence of matrix algebra renders the perturbative calculation of (TiP^^ and quantities built from 
it straightforward.) Thus, to sufficient accuracy, we obtain 

-i//g2|e-?V(ri+r2) r dsids2ds3dS4 e-?'(^3+.4)/2 ^-^in+r.)-' ^^' = 1 



2' 

^g^C-Ti+ra)"^ ^^^^ Tp{S3.p+S4.p 
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FIG. 1. Dependence of the gel fraction q on fi^: stable solution (full line); unstable solution (dotted line). 



57 



FIG. 2. Scaling function -iT{d) (full line) for the probability distribution of localisation lengths; asymptotic form for — > 
(dotted line); asymptotic form for ^ — » oo (broken line). 



FIG. 3. Scaling function u!{k) for the order parameter (full line); asymptotic form for — » (dotted line); asymptotic form 
for — » oo (broken line). 



TABLE I. Values of the order parameter in the liquid and amorphous solid states. 



Sector 


Order parameter 


Liquid state 
< 1) 


Amorphous solid state 
(/.^ > 1) 


one- replica 


(<3k)n + l 
{Qk)n+1 
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higher-replica 
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